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The  wind  bloweth  where  it  listeth 


and  thou  hearest  the  sound  thereof, 
but  canst  not  tell  whence  it  cometh 


and  whither  it  goeth 
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ABSTRACT 

V 

-  4 

Conventional  methods  for  monitoring  winds  aloft  have  inherent  spa¬ 
tial  and/or  temporal  limitations  that  severely  restrict  their  useful¬ 
ness.  An  alternate  solution  is  the  use  of  transhorizon  bistatic  radio 
techniques.  By  using  a  phase-coherent  system  with  some  form  of  beam- 
swinging  capability,  it  is  possible  to  measure  the  wind-produced  Doppler 
shifts  associated  with  moving  off-path  scatterers. 

In  this  experiment  a  rapid  data-gathering  system  is  used  to  re¬ 
motely  probe  the  wind-field  velocity  structure  in  the  common  volume  re¬ 
gion  of  a  164-km  transhorizon  propagation  path.  The  phase-stable  sys¬ 
tem  operates  at  3.2  GHz  and  employs  a  12-element  horizontal  receiving 
antenna  array.  The  amplitudes  and  relative  phases  of  the  output  sig¬ 
nals  from  the  individual  elements  are  sequentially  sampled  and  recorded 

!  I 

1  at  a  40-Hz  rate.  Off-line  processing  of  the  amplitude  and  phase  data  /  /*  .■ 
samples  is  used  to  '’scan”  in  azimuth  a  0.3°  wide  fan  beam  through  a 
3.88°  sector  and  to  obtain  virtual  simultaneous  Doppler  spectra  measure¬ 
ments  at  a  continuum  of  angles  within  the  sector  scanned.  An  algorithm 
developed  to  identify  the  apparent  source  direction  of  each  spectral 
component  permits  an  appropriate  weighting  to  be  applied  to  the  spectra 
on  the  basis  of  angl e-of-arrival .  The  average  Doppler  shift  (spectrum 
centroid)  as  a  function  of  the  array  pointing  angle  is  interpreted  in 
terms  of  the  average  transverse  wind  speed  and  direction.  Doppler  re¬ 
lations  applicable  to  previous  experiments  are  modified  fo\account  for 
the  influence  of  differences  in  the  longitudinal  path  geometry.  An 
analytical  expression  employing  a  path  asymmetry  factor  is  developed 
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for  the  wind-produced  Doppler  from  a  single  scatterer  model.  Wind 
velocities  inferred  from  the  average  Doppler  shift  data  are  compared 
with  in  situ  wind  observation  data  obtained  by  optically  tracking  pi 
lot  balloons  released  into  the  common  volume. 

.  Measurement  results  show  the  average  Doppler  shifts  to  have  an 
almost  linear  dependence  on  angl e-of-arrival .  Wind-produced  Doppler 
shifts  resulting  from  average  transverse  velocities  up  to  about  5  m/ 
have  been  measured  and  found  consistent  both  in  magnitude  and  direc¬ 
tion  with  the  prevailing  winds  observed  in  the  lower  strata  of  the 
common  volume  region.  On  numerous  occasions  Doppler  measurement  re¬ 
sults  reflected  much  higher  average  velocities  (-40  m/s)  than  could 
be  attributed  to  normal  low-altitude  winds.  Time-spatial  analyses 
of  these  anomalous  spectra  are  indicative  of  aircraft-induced  atmos¬ 
pheric  disturbances.  The  spectral  signatures  of  these  disturbances 
may  provide  additional  insight  into  the  generation  and  decay  of  air¬ 
craft  trailing  vortices. 
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Chapter  1 


INTRODUCTION 

This  dissertation  describes  the  results  of  a  microwave  transhorizon 
experiment  in  which  a  rapid  data-gathering  antenna  array  is  used  as  a 
remote  sensing  probe  of  low-altitude  tropospheric  winds. 

1 .1  Background 

The  initial  transhorizon  propagation  of  an  ultrahigh  frequency  sig¬ 
nal  was  reported  by  Karconi  in  1932  after  he  successfully  transmitted  a 
500-MHz  signal  over  a  distance  of  168  miles  (270  km).  The  underlying 
physical  principles  involved  in  transhorizon  propagation  subsequently 
became  a  subject  of  discussion  and  some  controversy  among  researchers. 
The  smooth-earth  diffraction  theory  (Watson,  1918  and  1919;  Van  der  Pol 
and  Bremmer,  1937,  1938,  and  1939),  which  suggests  an  exponential  de¬ 
crease  in  the  received  power  with  distance,  was  the  generally  accepted 
explanation  for  transhorizon  propagation  in  the  early  1940s.  With  the 
advent  of  high-power  transmitters  and  sensitive  receivers  during  the 
war,  signal  reception  far  beyond  the  horizon  became  more  frequent  (Bul- 
lington,  1955).  The  received  signal  levels,  however,  were  consistently 
found  to  be  far  above  the  levels  predicted  by  the  diffraction  theories. 
Concerted  efforts  were  made  by  investigators  in  the  late  1940s  to  find 
a  plausible  explanation  for  the  transhorizon  signal  observations.  Pe- 
keris  (1947)  suggested  the  signals,  based  on  their  appearance,  were  a 
result  of  some  scattering  process  from  refractive- index  variations  in 
the  atmosphere.  As  a  consequence,  several  theoretical  propagation 


1 


models  evolved  that  took  into  account  the  prevailing  atmospheric  struc¬ 
ture.  The  preponderance  of  experimental  data  taken  over  the  years  sup¬ 
ports  the  validity  of  two  basic  theoretical  models,  one  derived  from  a 
turbulence  scattering  theory  (Booker  and  Gordon,  1950;  Staras,  1952) 
and  the  other  derived  from  a  layer  reflection  theory  (Bauer,  1956;  Friis 
et  al.,  1957): 

(1)  In  the  turbulence  theory  the  propagation  phenomena 
is  characterized  as  a  scattering  process  resulting 
from  turbulent  conditions  in  the  atmosphere  that 
cause  random  irregularities  in  the  refractive  pro¬ 
perties  of  the  medium. 

(2)  In  the  layer  theory  the  propagation  phenomena  is 
characterized  as  a  reflection  process  resulting 
from  stratified  layers  in  the  atmosphere  that  pro¬ 
duce  sharp  discontinuities  in  the  refractive  index 
gradient. 

Much  of  the  experimental  and  theoretical  work  conducted  since  the 
introduction  of  the  two  fundamental  theories  has  dealt  with  the  study 
of  the  physical  mechanisms  supporting  propagation  beyond  the  horizon. 

The  early  experimental  research,  however,  was  influenced  to  a  large  ex¬ 
tent  by  the  demands  for  reliable  long-distance  radio  communications 
(Bullington,  1950;  Mellen  et  al.,  1955).  Consequently,  one  of  the 
principal  concerns  of  researchers  was  the  signal  fading  problems  often 
experienced  on  transhorizon  paths.  The  cumulative  results  of  the  nu¬ 
merous  experiments  conducted  up  to  about  1960  provided  sufficient  data 
to  establish  a  statistical  basis  for  predicting  (Baris  et  al.,  1962) 
the  average  expected  level  of  performance  of  transhorizon  radio  trans¬ 
mission  systems. 

By  1960  the  tropospheric  propagation  characteristics  beyond  the 


horizon  were  fairly  well  understood.  Investigators  subsequently  began 
to  look  into  ways  of  distinguishing  between  the  propagation  mechanisms 
associated  with  the  turbulence  and  layer  model  theories.  At  Stanford 
researchers  (Waterman  et  al.,  1957;  Waterman,  1958;  Lee,  1961;  Stroh- 
behn,  1963;  Waterman  and  Strohbehn,  1963;  Strohbehn  and  Waterman,  1966; 
Cox,  1967;  Cox  and  Waterman,  1969  and  1971;  Cianos,  1971;  Cianos  and 
Waterman,  1973;  Waterman,  1973;  Cianos,  1978)  have  used  various  experi¬ 
mental  techniques  to  differentiate  between  the  two  prevailing  atmos¬ 
pheric  structures  suggested  by  the  theoretical  models.  These  empirical 
methods  provided  in  effect  a  means  of  remotely  probing  the  atmosphere 
with  radio  waves  —  a  technique  in  which  analyses  of  specific  parameters 
of  the  received  signal  are  used  to  infer  certain  characteristics  about 
the  intervening  medium. 

Several  important  meteorological  parameters  can  be  measured  re¬ 
motely  using  transhorizon  radio  techniques.  Among  the  many  measureable 
parameters,  the  wind  velocity  structure  in  the  lower  atmosphere  has  re¬ 
ceived  considerable  attention  by  investigators  (Atlas  et  al.,  1968; 
Birkemeier  et  al . ,  1968;  Atlas  et  al . ,  1969;  Atlas,  1969;  Birkemeier  et 
al.,  1969;  Lammers  and  Olsen,  1973;  Chadwick  et  al . ,  1975)  over  the 
past  decade.  Larly  investigations  of  the  signal  fading  phenomena  on 
transhorizon  paths  have  demonstrated  marked  correlations  between  signal 
fade  rates  and  local  transverse  wind  conditions  (Laaspere,  1958;  Craw¬ 
ford  et  al . ,  1959;  Doherty,  1959;  Gjessing,  1964).  Refractive  struc¬ 
tures  within  the  common  volume  drift  with  the  mean  wind  velocity  which 
results  in  stochastic  phase  variations  in  the  received  signal.  These 
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phase  variations  are  reflected  in  the  Doppler  spectrum  and  are  depen¬ 
dent  upon  the  spatial  distribution,  size,  and  relative  velocities  of 
the  moving  atmospheric  structures.  The  analytical  and  experimental 
work  performed  in  this  area  have  provided  the  motivation  for  the  ex¬ 
periment  described  in  this  dissertation  which  seeks  to  further  our  un¬ 
derstanding  of  the  spatial  and  temporal  variations  of  atmospheric  mo¬ 
tion  produced  by  the  wind. 

1 .2  Significance  of  Experiment 

Knowledge  of  wind  information  in  localized  regions  of  the  lower 
atmosphere  is  of  vital  importance  not  only  to  meteorologists  but  to 
researchers  in  such  diverse  areas  as  air- traffic  safety,  environmental 
pollution,  energy,  and  even  agriculture.  The  conventional  rawinsonde 
technique  of  monitoring  winds  aloft  has  inherent  spatial  and  temporal 
limitations  that  severely  restrict  its  usefulness.  Situations  requir¬ 
ing  near  real-time  wind  information  in  selected  regions  of  the  atmo¬ 
sphere  are  not  adequately  served  by  the  widely  spaced  and  infrequent 
(every  12  hours)  launches  of  rawinsondes.  A  reduction  in  either  the 
time  interval  or  the  spacing  between  launches  becomes  economically  pro¬ 
hibitive  in  terms  of  the  increased  manpower  and  instrumentation  re¬ 
quirements. 

An  alternate  means  of  wind  sounding  in  the  lower  atmosphere  has 
recently  been  proposed  by  Chadwick  et  al.  (1975).  The  remote  probing 
technique  proposed  involves  the  use  of  a  network  of  bistatic  radio 
transmission  systems.  The  proposal  stems  from  the  results  of  experi¬ 
mental  work  by  Birkemeier  et  al.  (1968  and  1969)  and  Laimners  and  Olsen 
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(1973).  Analyses  by  Atlas  (1969)  and  Atlas  et  al.  (1965,  1968,  and 
1969)  have  provided  considerable  insight  into  the  elementary  theory  and 
have  enabled  the  experimental  observations  to  be  interpreted  in  meteoro¬ 
logically  significant  terms.  Observed  systematic  increases  in  the  aver¬ 
age  Doppler  shift  with  offset  angle  from  the  great-circle  bearing  have 
been  attributed  to  the  wind-produced  drift  of  refractive  structures  in 
the  atmosphere.  By  incorporating  narrow-beam  antennas  with  some  form  of 
beam-swinging  capability,  it  is  possible  to  remotely  monitor  the  trans¬ 
verse  (cross-path)  component  of  the  wind  over  a  wide  spatial  extent  using 
a  bistatic  phase-coherent  transhorizon  system.  Aligning  two  such  systems 
orthogonal  to  one  another  permits  a  complete  characterization  of  the  hori 
zontal  wind-field. 

Previous  experiments  have  employed  synchronously  offset  steerable 
antennas  to  probe  selected  regions  in  the  atmosphere  and  have  been  con¬ 
ducted  on  longitudinally  symmetrical  propagation  paths.  These  are  trans¬ 
horizon  paths  where  the  common  volume  region  is  symmetrical  about  the  mid 
path  plane  between  the  receiver  and  transmitter  terminals.  Over  highly 
asymmetrical  paths  the  common  volume  region  is  skewed  significantly  to 
one  side  of  the  mid-path  plane.  The  resultant  Doppler  relations  must 
then  reflect  the  influence  of  the  longitudinal  position  of  the  dominant 
scattering  region.  Analytical  results  applicable  to  both  longitudinally 
symmetrical  and  asymmetrical  paths  are  presented  in  this  dissertation 
along  with  experimental  measurement  data.  Unique  features  of  the  experi¬ 
ment  pertinent  to  the  remote  probing  of  tropospheric  winds  are 

(1)  Angl e-of-arrival  measurements  of  wind-produced  Doppler 
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are  derived  from  processing  array-sampled  amplitude 
and  relative  phase  data. 

(2)  Rapid  sampling  features  of  the  data-gathering  array 
permit  the  simultaneous  measurement  of  the  upwind 
and  downwind  Doppler  spectra. 

(3)  Measurements  are  made  on  a  longitudinally  asymmetrical 
propagation  path  where  the  principal  scattering  region 
is  centrally  located  between  the  receiver  and  the  mid¬ 
path  plane. 

1 .3  Description  of  Experiment 

To  investigate  the  spatial  and  temporal  variations  of  tropospheric 
winds,  a  42-hour  transhorizon  propagation  experiment  was  conducted  over 
a  164-km  asymmetrical  path  at  3.2  GHz.  The  propagation  path  extended 
from  Jackson  Butte,  California  (transmitter)  to  the  foothills  behind 
Stanford  University  (receiver).  A  rapid  data-gathering  system  (Cox  and 
Waterman,  1969)  used  in  previous  transhorizon  experiments  (Cox,  1967; 

Cox  and  Waterman,  1971;  Cianos,  1971;  Cianos  and  Waterman,  1973;  Water¬ 
man,  1973;  Cianos,  1978)  was  modified  to  remotely  probe  the  wind-field 
structure  in  the  common  volume  region  above  the  Livermore  Valley  about 
40  km  from  the  receiver. 

In  this  experiment  the  phase-coherent  system  employs  a  horizontal 
12-element  linear  receiving  antenna  array.  The  amplitudes  and  relative 
phases  of  individual  elements  are  sequentially  sampled  and  recorded  at 
a  40-Hz  rate.  An  effective  scan  period  of  5  ms  is  used  to  acquire  the 
data  from  all  twelve  elements.  Off-line  signal  processing  is  performed 
on  the  amplitude  and  phase  data  samples  to  "scan"  in  azimuth  a  0.3°  wide 
beam  through  a  3.88°  sector  and  to  obtain  corresponding  angle-of-arrival 
Doppler  spectra.  The  average  Doppler  shift  (spectrum  centroid)  is 
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computed  as  a  function  of  the  array  pointing  angle  and  then  interpreted 
in  terms  of  the  average  transverse  wind  speed  and  direction.  Results 
are  compared  with  in-situ  wind  observation  data  obtained  by  optically 
tracking  pilot  balloons  released  into  the  common  volume. 

The  experimental  data  presented  in  this  dissertation  were  obtained 
from  tests  conducted  on  April  21-22,  1978.  Ten  minutes  of  transhorizon 
radio  data  were  recorded  each  hour.  A  pilot  balloon  was  launched  hour¬ 
ly  between  0800  and  1800  hours  in  the  middle  of  every  10-minute  data 
recording  period. 


1 .4  Summary  of  Results 

Some  of  the  significant  results  of  this  research  are: 

(1)  Wind-Ooppler  measurement  results  compare  most  favorably 
with  in-situ  PIBAL  data  in  lower  strata  of  the  common 
volume. 

(2)  Spectral  broadening  of  wind-produced  spectra  are  ob¬ 
served  during  afternoon  hours  when  the  atmosphere  is 
well  mixed. 

(3)  At  symmetrical  off-axis  angles,  measured  upwind  and 
downwind  Doppler  spectra  are  observed  to  be  nearly 
antisymmetric  about  the  zero-Doppler  frequency.  This 
suggests  some  degree  of  spatial  uniformity  in  the 
transverse  winds. 

(4)  On-axis  the  Doppler  spectra  are  generally  symmetrical 
about  zero-Doppler  suggesting  the  absence  of  vertical 
wind  components  and/or  spatial  non-uniformities  in  the 
transverse  wind  velocity. 

(5)  Aircraft  trailing  vorticies  produce  spectral  broaden¬ 
ing  and  spatial  asymmetries  in  the  ambient  wind  pro¬ 
duced  spectra. 

(6)  Spectral  and  transit-time  analyses  of  aircraft-induced 
Doppler  permit  a  determination  of  the  spatial  extent, 
range,  and  velocity  of  trailing  vorticies. 
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(7)  Analytical  results  show  the  Doppler  shift  for  a 
scatterer  at  a  fixed  cross-path  position  to  be 
greater  on  a  longitudinally  asymmetrical  path  vis 
a  vis  a  symmetrical  path  of  equal  path  length. 
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Chapter  2 


THEORETICAL  DISCUSSION 

Signals  received  via  transhorizon  tropospheric  scattering  proces¬ 
ses  are  frequently  subject  to  rapid  fading  and  spectral  instabilities. 
This  is  a  consequence  of  the  interference  nature  of  the  Doppler-shifted 
components  associated  with  moving  off-path  scatterers.  Interference 
patterns  are  set  up  at  the  receiver  by  signal  components  scattered  off 
refractive  eddy  structures  transported  by  the  wind  through  regions  of 
the  atmosphere  mutually  illuminated  by  the  transmitting  and  receiving 
antenna  beams.  As  suggested  by  Atlas  et  al .  (1969),  atmospheric  struc¬ 
tures  on  the  upwind  side  of  the  beam  produce  positive  Doppler  frequency 
shifts  that  beat  with  negative-shifted  Doppler  components  associated 
with  eddies  on  the  downwind  side  of  the  beam.  These  Doppler  shifts 
progressively  increase  with  the  off-path  displacement  of  the  scattering 
region  and  have  been  shown  (Birkemeier  et  al.,  1968;  Atlas  et  al.,  1969) 
to  be  directly  proportional  to  the  cross-path  wind  speed  in  the  lower 
segments  of  the  common  volume.  A  single-scatterer  model  (Birkemeier 
et  al.,  1968)  can  be  used  to  illustrate  the  physical  concepts  involved 
and  to  show  the  influence  of  the  path  geometry  and  source  position  on 
the  basic  Doppler  relations. 

2.1  Single  Scatterer  Model 

In  the  analysis  presented  in  this  section  a  single  scatterer  is 
assumed  to  be  responsible  for  the  signal  observed  at  the  receiver.  As 
shown  in  Figure  2.1,  the  point  source  scatterer  is  located  at  the 


intersection  of  two  perfectly  narrow  antenna  beams.  The  Cartesian  co¬ 
ordinate  system  used  In  this  model  has  its  origin  centrally  located 

l 

along  the  chordal  axis  between  the  two  radio  terminals.  The  x,  y,  and 

z  coordinates  define  respectively  the  longitudinal,  transverse,  and 

vertical  spatial  positions  of  the  scatterer.  The  total  ray-path  length, 

L  ,  from  the  transmitter  to  the  receiver  is  the  sum  of  the  respective 

path  lengths,  L.  and  L  ,  from  the  transmitter  to  the  scatterer  and 
t  r 

from  the  receiver  to  the  scatterer.  Denoting  the  chordal  distance  be¬ 
tween  the  terminals  as  2D  ,  the  total  ray-path  length  can  be  expres¬ 
sed  in  terms  of  the  position  coordinates  of  the  scatterer  as 

L  =  C  ( D-t-x ) 2  +  y2  +  z2]*5  +  [(D-x)2  +  y2  +  z2]*5  (2.1 ) 


If  we  assume  a  4/3  effective  earth  radius  (Bean  and  Dutton,  1966)  to  ac¬ 
count  for  the  refractive  bending  of  radio  rays  in  the  lower  atmosphere, 
we  may  use  straight  lines  for  the  ray  paths  in  applying  the  above  expres¬ 
sion. 

It  is  well  known  that  a  moving  source  produces  a  Doppler  shift, 
fd  ,  at  the  receiver  that  is  proportional  to  the  time  rate  of  change  of 
the  ray-path  length.  This  Doppler  principle  can  be  expressed  as 


f  =  .L4k 

Td  X  dt 


(2.2) 


where  X  is  the  wavelength  of  the  transmitted  signal.  It  is  seen  that 
if  the  effective  path  length  decreases  (increases)  with  time,  a  positive 
(negative)  Doppler  shift  will  be  observed  at  the  receiver.  The  Doppler 
shift  can  also  be  interpreted  as  the  rate  at  which  the  scatterer  moves 


across  surfaces  of  constant  path  length.  Contours  of  constant  path 
length  are  ellipsoidal  surfaces  (Birkemeier  et  al.,  1968)  with  the  radio 
terminals  as  foci.  If  the  scatterer  moves  with  a  velocity,  7$  ,  an 

equivalent  Doppler  expression  can  be  written  as 


f 


d 


lv 

X  vs 


•  VL 


(2.3) 


Upon  substituting  the  expression  given  in  (2.1)  for  the  total  ray-path 
length,  the  Doppler  shift  from  a  single  moving  scatterer  becomes  expres¬ 
sed  as 


[ux  +  vy  +  wz]L  +  uD[(y2-2xD)^  -  (y2+2xP)^] 


(y2+2xD)is  (y2-2xD)Js 


(2.4) 


where  In  the  above  y2  =  [x2  +  y2  +  z2  +  D2]  and  u,  v,  and  w  are  re¬ 
spectively  the  longitudinal,  transverse,  and  vertical  velocity  components 
of  the  moving  scatterer.  From  the  general  result  in  (2.4)  specific  Dop¬ 
pler  expressions  for  symmetrical  and  asymmetrical  paths  will  be  derived 
to  illustrate  the  influence  of  the  scatterer's  longitudinal  position. 


2.2  Symmetrical  Path  Case 

On  longitudinally  symmetrical  paths,  the  scattering  regions  pri¬ 
marily  responsible  for  the  received  signal  are  located  in  close  proxi¬ 
mity  of  the  mid-path  plane.  As  a  result,  a  moving  point-source  scat¬ 
terer  has  relatively  small  longitudinal  displacements,  i.e.  x  «  D  and 
is  generally  confined  to  scattering  regions  in  the  common  volume  where 
D2  »  [x2  +  y2  +  z2]  .  Applying  a  first-order  approximation  to  the  terms 
(y2+2xD)15  and  (y2-2xD)^  in  (2.4),  the  Doppler  shift  produced  by  a 
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moving  scatterer  on  a  symmetrical  path  becomes  expressible  as 

(2.5) 

This  result  arrived  at  by  Birkemeier  et  al .  (1968)  has  been  applicable 
to  previously  reported  transhorizon  wind-sensing  experiments  since  such 
experiments  have  heretofore  been  conducted  on  near-symmetrical  paths. 

Equation  (2.5)  can  be  further  simplified  by  considering  the  geome¬ 
try  of  the  symmetrical  path  profile  in  Figure  2.2  and  the  relative  con¬ 
tributions  of  the  wind  velocity  components.  It  can  be  shown  that  the 
term  (1  -4D2/L2 )  is  small  and  equivalent  to  4z2/L2  =  4h2/D2  where  zm 
and  hm  are  respectively  the  mid-path  heights  of  the  scatterer  above  the 
chordal  axis  and  the  earth's  surface.  Since  the  longitudinal  displace¬ 
ment,  x  ,  of  the  scatterer  is  also  small,  the  Doppler  contribution  from 
the  longitudinal  wind  can  be  neglected,  the  longitudinal  motion  produ¬ 
ces  negligible  Doppler  effects  because  of  the  elongated  shapes  of  the 
ellipsoidal  surfaces  of  constant  path  length.  A  scatterer  moving  in 
the  longitudinal  direction  cuts  across  the  surfaces  at  a  much  slower 
rate  than  one  moving  in  the  transverse  or  vertical  directions.  The  ef¬ 
fect  of  the  vertical  wind  can  in  most  instances  be  neglected  since  its 
magnitude  is  normally  much  less  than  that  of  the  horizontal  wind.  It 
then  becomes  evident  that  the  Doppler  shift  on  a  symmetrical  path  is 
linearly  proportional  to  the  cross-path  displacement,  y  ,  of  the  mov¬ 
ing  scatterer  and  may  be  explicitly  given  as 

Td  XL  (2.6) 


4  ux(1-4D2/L2)  +vy  +  wz 


1  -  1  6x2D2/Li‘ 
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2.2.  SYMMETRICAL  PATH  PROFILE. 
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The  Doppler  shift  may  also  be  expressed  in  terms  of  the  azimuthal  angu¬ 
lar  position,  a  ,  of  the  scatterer.  The  Doppler  relationship  is  then 
written  as 

fd*"rsin(ct)  (2'7) 


which,  for  small  angles,  can  be  approximated  with  negligible  error  as 


2v« 

X 


(2.8) 


2.3  Asymmetrical  Path  Case 

Consider  now  the  path  profile  shown  in  Figure  2.3  where  the  prin¬ 
cipal  scattering  region  in  the  common  volume  is  skewed  significantly  to 
one  side  of  the  mid-path  plane.  We  begin  this  treatise  by  again  looking 
at  the  terms  (y2+2xD)^  and  (y2-2xD)^  in  the  general  Doppler  shift  ex¬ 
pression  given  in  (2.4)  but  view  them  from  the  standpoint  of  an  asym¬ 
metrical  path.  For  this  situation,  it  is  assumed  that  the  longitudinal 
displacement,  x  ,  is  such  that  (D-x)2  »  y2  +  z2  for  all  values  of 
x  ,  y  ,  and  z  within  the  common  volume.  With  this  constraint,  a  first- 
order  approximation  to  the  terms,  (y2+2xD)^  and  (y2-2xD) 2  ,  results  in 
a  Doppler  shift  on  an  asymmetrical  path  that  can  be  expressed  as 


L 

XD2 


2ux(l -2D/L)  +  vy  +  wz 


1  -  x2/D2 


(2.9) 


This  result  reduces  to  that  for  the  symmetrical  path  case  when  x  «  D  . 
Hence,  the  Doppler  shift  expression  given  in  (2.9)  is  a  more  general  re¬ 
sult  that  is  applicable  to  both  symmetric  and  asymmetric  path  geometries. 
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A  simplification  of  (2.9)  can  be  made  by  introducing  a  path  asym¬ 
metry  factor  resulting  from  the  terrain  profile  along  the  path  and  the 
relative  heights  of  the  radio  terminals.  The  asymmetry  factor,  s  ,  is 
defined  as  the  ratio  of  the  longitudinal  displacement  of  the  scattering 
source  to  the  chordal  distance  from  the  mid-path  plane  to  either  of  the 
radio  terminals,  i.e.  s  =  x/D  .  It  is  zero  for  a  symmetrical  path. 

When  the  scattering  region  is  skewed  to  one  side  of  the  mid-path  plane, 
however,  the  magnitude  of  the  asymmetry  factor  becomes  greater  than  zero 
but  always  remains  less  than  unity.  Using  this  factor,  the  Doppler  shift 
becomes 


L  us(°Vat  +  var  +  wer 

xd  :  77  +  T~ : 


1  -  ss 


1  +  s 


(2.11) 


where  small  angle  approximations  with  respect  to  the  receiver  and  trans 
mltter  are  used,  i.e.  ar  =  ,  ot  «  5^  ,  , 

and  .  For  equal  transverse  and  longitudinal  wind  veloci¬ 

ties,  the  longitudinal  Doppler  contribution  again  remains  relatively 
small  and  can  be  neglected.  If  the  vertical  wind  contribution  is  also 


neglected  and  the  first-order  approximation,  L  =  2D  ,  is  substituted 
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for  the  ray-path  length,  the  Doppler  shift  expression  (2.10)  simplifies 
to 


which  can  be  rewritten  in  terms  of  the  azimuthal  angular  offset,  ar  , 
from  the  receiver  as 


2va 

f d  “  '  XlT+iT 


(2.13) 


Comparing  the  above  results  with  those  for  the  symmetrical  path  case, 
equations  (2.6)  and  (2.8),  it  is  seen  that  the  Doppler  shift  for  a  scat- 
terer  at  a  given  cross-path  displacement  is  greater  on  an  asymmetrical 
path  than  on  a  symmetrical  path.  This  is  consistent  with  the  interpreta¬ 
tion  of  the  Doppler  shift  as  the  rate  at  which  a  scatterer  cuts  across 
surfaces  of  constant  path  length.  The  separation  between  these  ellipsoi¬ 
dal  surfaces  is  greatest  in  the  mid-path  plane  and  becomes  smaller  to¬ 
wards  either  of  the  two  terminals.  Consequently,  a  scatterer  moving 
through  the  common  volume  region  on  an  asymmetrical  path  cuts  across  the 
equiphase  surfaces  at  a  more  rapid  rate  than  on  a  symmetrical  path. 

On  the  asymmetrical  Jackson-Stanford  path  used  in  the  experiment 
discussed  in  this  dissertation,  the  principal  scattering  region  is  in  the 
vicinity  of  a  transverse  plane  about  40  kilometers  from  the  receiver. 

The  average  longitudinal  displacement,  x  ,  of  a  scatterer  on  this  path 
is  such  that  s  =  H  •  The  analytical  expression  for  the  Doppler  shift 
applicable  to  this  path  is 
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(2.14) 


where  a  Is  measured  in  radians.  If 
r 

ted  to  degrees  and  the  numerical  value 
\  ,  a  more  simplified  expression  for 
Stanford  path  can  be  written  as 


va 


f  J  *  - 


the  angle  measurement  is  conver- 
(9.4  x  10"2m)  is  substituted  for 
the  Doppler  shift  on  the  Jackson- 


(2,15) 


where  in  the  above  expression  v  is  units  of  m/s. 

2.4  Multiple  Scattering  Effects 

In  the  previous  sections,  a  single  scattering  source  was  assumed  and 
as  a  result,  Doppler  effects  were  limited  to  single  frequency  shifts.  A 
multiplicity  of  scatterers  randomly  distributed  throughout  an  illuminated 
region  in  the  atmosphere  produces  a  spectrum  of  Doppler  shifts.  The 
spectral  composition  of  the  Doppler  shifts  is  dependent  upon  the  spatial 
distribution  and  the  relative  velocities  of  the  individual  scatterers. 
Each  scatterer  gives  rise  to  a  discrete  Doppler  shift  in  accordance  with 
the  analysis  developed  earlier.  A  number  of  the  scatterers,  however,  may 
contribute  to  the  same  spectral  component.  The  extent  of  the  contribu¬ 
tions  from  the  many  off-path  scatterers  is  reflected  in  the  shape  of  the 
Doppler  spectrum,  which  is  dependent  upon  the  relative  strength  of  the 
received  signal  as  a  function  of  the  scatterers'  cross-path  positions. 

The  signal  power  received  from  a  particular  azimuthal  direction  is  di¬ 
rectly  proportional  to  the  product  of  the  net  antenna  gain  and  the  re¬ 
flectivity  or  scatter  cross-section  of  the  medium.  Thus,  the  Doppler 
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spectrum  arising  from  a  distribution  of  scatterers  is  dependent  upon  the 
antenna  beam  configuration  and  the  basic  propagation  mechanism. 

An  average  Doppler  shift,  f^  ,  can  be  attributed  to  a  mean  flow 
or  drift  of  the  atmospheric  scatterers.  Birkemeier  et  al .  (1968)  have 
specified  analytically  the  average  Doppler  shift  for  antennas  synchro¬ 
nously  aligned  in  azimuth  at  an  angle,  a  ,  as 

d 

7d  “  -  r/“G(“  -aa)W(a)da  //g(  a-aa)W(a)da  (2.16) 

where  G(a-a  )  is  the  combined  antenna  gain  function  and  W(ot)  is  the 

d 

scattered  power  per  unit  azimuthal  angle,  this  result,  however,  applies 
to  a  symmetrical  path.  A  more  general  result  incorporates  the  path  asym¬ 
metry  factor,  s  ,  and  can  be  given  in  terms  of  the  angular  offset, 

a  ,  measured  from  the  receiver  location,  as 
r 

?d  =  yaG(a-ar)W(a)da  /J G(a-ar )W(a)da  (2.17) 

In  these  expressions  for  the  average  Doppler  shift,  various  atmospheric 
models  can  be  applied  to  account  for  the  angular-dependence  scattering 
function,  W(a)  .  The  classical  isotropic  turbulence  model  (Tatarski, 
1961)  predicts  a  minus  11/3  power-law  variation  in  the  scattering  cross- 
section  as  a  function  of  scattering  angle,  9  .  The  scattering  angle  is 

defined  as  the  angle  between  the  transmitted  and  received  ray-path  di¬ 
rections.  For  the  isotropic  turbulence  model,  the  scattering  cross- 

_  li 

section  is  proportional  to  sin  3  (9/2)  .  Birkemeier  et  al .  (1977) 

have  shown  that  for  small  scattering  angles  near  the  mid-path  plane 
9  =  Z^a2  +  32  .  This  expression  for  the  scattering  angle  is  valid  for 
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near-symmetric  propagation  paths.  However,  a  more  general  result  that 
Is  applicable  to  asymmetrical  paths  as  well  shows  the  scattering  angle 
to  be  given  approximately  as 

®  *  >/<*£  +  ef  +  +  (2,18) 

Using  the  small  angle  approximations  developed  earlier,  the  above  expres¬ 
sion  can  be  written  as  a  function  of  the  path  asymmetry  factor,  s  ,  in 
terms  of  the  angular  displacements,  and  3^  ,  from  the  receiver  as 

6  “  t4t  V«*r  +  ej  (2.19) 

An  equivalent  expression  for  the  Jackson-Stanford  path,  which  has  a  mini¬ 
mum  scattering  angle  of  0.84°  occurring  at  an  elevation  angle  of  0.25°, 
is 

9  «  0.5  +  |  +  B J  (2.20) 


where  6  >  0.25°  . 

r  — 

In  addition  to  the  azimuthal  variations,  the  above  relationships 
show  the  extent  of  the  influence  of  the  vertical  positions  of  the  scat- 
terers  on  the  scattering  function.  The  average  Doppler  shift  from  a 
given  azimuthal  direction  is  a  spatially  averaged  result  of  the  trans¬ 
verse  motions  of  scatterers  at  various  altitudes  within  the  common 
volume.  The  wind-produced  drift  of  the  scatterers  generally  increases 
with  altitude.  However,  the  turbulence  model  assumed  produces  a  rapid 
decrease  in  the  scattered  power  with  elevation  angle.  This  suggests 
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the  average  Doppler  shift  is  primarily  associated  with  the  motion  of 
scatterers  in  the  lower  segments  of  the  illuminated  scattering  region. 
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Chapter  3 


DESCRIPTION  OF  EXPERIMENT 


3.1  Propagation  Path 

The  experiment  described  in  this  report  was  conducted  over  the  pro¬ 
pagation  path  shown  in  Figure  3.1.  The  164-km  (102  miles)  path  is  align¬ 
ed  along  a  SW-NE  (230°  -  50°)  azimuthal  bearing  between  Jackson  Butte 
(transmitter  site)  and  Stanford,  California  (receiver  site).  The  trans¬ 
mitter  was  situated  on  a  725-meter  (2375  ft)  peak  just  east  of  Jackson, 
California  and  the  receiver  was  located  110  meters  (360  ft)  above  sea  le¬ 
vel  in  the  foothills  behind  Stanford  University.  A  topographical  profile 


Fig.  3.1.  PROPAGATION  PATH  AND  SURROUNDING  AREA. 
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Fig.  3.2.  TOPOGRAPHICAL  PROFILE  OF  PROPAGATION  PATH. 


(Figure  3.2)  of  the  path  shows  the  existence  of  three  mountain  ridges 
that  are  nearly  perpendicular  to  the  path  between  the  two  terminals. 
These  ridges  determine  the  local  radio  horizons  of  the  transmitting  and 
receiving  antenna  beams  which  are  aligned  so  as  to  mutually  illuminate  a 
region  from  about  425  meters  (1400  ft)  to  approximately  2  km  (6600  ft) 
above  sea  level.  The  path  geometry  introduces  a  minimum  scattering  an¬ 
gle  of  0.84° 

3.2  Instrumentation 

3.2.1  Data-Gathering  Array 

Angle-of-arrival  Doppler  measurements  are  derived  from  the  indivi¬ 
dual  amplitude  and  relative  phase  sampled  outputs  of  a  12-element  data 
gathering  antenna  array.  The  array,  which  operates  at  3.2  GHz,  is  shown 
in  Figure  3.3  and,  except  for  the  horizontal  orientation  and  polariza¬ 
tion,  is  identical  to  the  vertical  array  used  in  previous  experiments  at 
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Stanford  (Cox,  1967;  Cox  and  Waterman,  1969  and  1971;  Cianos,  1971; 

Clanos  and  Waterman,  1973;  Waterman,  1973;  and  Cianos,  1978).  The  data- 
gatherlng  array  is  aligned  nearly  broadside  to  the  propagation  path  and 
is  structurally  supported  beneath  the  2nd,  6th,  and  10th  elements.  Each 
element  is  horizontally  polarized  and  has  a  5°  half-power  beamwidth  which 
determines  the  vertical  extent  of  the  array  fan  beam.  When  uniformly  il¬ 
luminated,  the  array  produces  in  azimuth  a  0.3°  wide  beam  which  can  be 
scanned  through  a  3.88°  azimuthal  sector.  A  summary  of  the  physical  and 
operational  characteristics  of  the  array  is  outlined  in  Table  3.1. 


Table  3.1 

SUMMARY  CHARACTERISTICS  OF  DATA-GATHERING  ARRAY 


Physical 

Characteristics 

Number  of  elements 

Element  spacing 

Element  diameter 

Effective  aperture  length 

Twel ve 

1.4  m  (14.8  wavelengths) 

1.2  m  (13.0  wavelengths) 

15.2  m  (162.5  wavelengths) 

Operating  frequency 

3.2  GHz 

Operational 

Polarization 

Horizontal 

Characteristics 

Half-power  beamwidth* 

0.3°  AZ/5°  EL 

Grating  lobe  separation* 

3.88° 

Sidelobe  levels* 

<  -  13  dB 

*  for  a  uniformly  illuminated  conventional  additive  array 
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3.2.2  Array-Receiver  Channels 

Twelve  separate  receiving  channels  connect  to  the  twelve  separate 
elements  in  the  array.  A  typical  receiving  channel  is  illustrated  in 
Figure  3.4.  Each  channel  receives  simultaneously  the  transhorizon  sig¬ 
nal  at  3.2001  GHz  and  a  stable  1 ine-of-sight  reference  signal  at  3.2000 
GHz.  A  3.2300-GHz  local  oscillator  signal  (common  to  all  12  channels) 
is  mixed  with  both  the  transhorizon  and  reference  signals  in  mixer- 
preamps  mounted  directly  behind  each  array  element.  The  down-converted 
transhorizon  (29.9-MHz)  and  reference  (30.0-MHz)  signals  are  separated 
by  2-kHz  wide  crystal  filters.  The  30-MHz  IF  reference  is  amplified 
and  then  used  as  a  second  local  oscillator  which  mixes  with  the  29.9-MHz 
transhorizon  IF  signal.  The  resultant  100-kHz  second  IF  signal  is  nar¬ 
row-band  filtered  with  a  23-Hz  wide  crystal  filter  which  effectively 
sets  the  noise  bandwidth  of  the  receiving  channel.  After  further  ampli¬ 
fication  the  narrow-band  100-kHz  signal  is  passed  through  a  variable 
phase  shifter  and  a  variable  gain  amplifier  which  are  adjusted  and  set 
during  calibration.  This  signal  is  then  continuously  available  from 
each  array-element  receiving  channel. 

Sampling  of  the  twelve  100-kHz  signals  is  performed  by  a  computer- 
controlled  electronic  commutator.  A  block  diagram  of  the  amplitude  and 
phase  sampling  equipment  is  shown  in  Figure  3.5.  Every  25  ms  the  com¬ 
mutator  is  sequentially  stepped  through  all  12-element  positions.  The 
commutator  takes  five  milliseconds  to  completely  scan  all  twelve  ele¬ 
ments.  The  commutator  output  is  bandpass  filtered  to  reduce  switching 
transients.  The  sample  100-kHz  output  is  phase  compared  with  a  stable 
100-kHz  reference  standard.  Additionally,  the  signal  is  applied  to  two 
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Fig.  3.4.  BLOCK  DIAGRAM  OF  AN  ARRAY  ELEMENT  RECEIVING  CHANNEL. 

Parameters 


1st  Mixers  (part  of  LEL  type 
SBC-6 

Noise  figure:  9  to  10  dB 
(depends  on  unit) 

30-MHz  IF  Preamplifiers  (part 
of  LEL  type  SBC-6) 

Gain:  25  to  30  dB  RF  to  IF 
(depends  on  unit) 

30-MHz  IF  Amplifiers  (LEL  type 

ITA-ir 

Adjustable  gain:  20  to  70  dB 
Bandwidth:  3  MHz 

29.9-MHz  and  30,0-MHz  Filters 
(Hughes  special  design) 

Bandwidth:  2  kHz 
Quartz  crystals 

2nd  Mixers 

Balanced  diodes  for  LO  (refer¬ 
ence  signal ) 

Noise  Suppression 


lOO-kHz  Filters  (Blackhawk  Networks 
Corporation) 

Bandwidth:  23  Hz 
Quartz  Crystals 

Phase  Shifters  (Nil sen  V42-01  Vari- 
ogon) 

Adjustable  0°  to  360°  continuously 
Loss  independent  of  phase  setting 

100-kHz  Output  Amplifiers 


Low  output  impedance:  5  to  10  n 
Phase  shift  independent  of  gain 
setting 

Gain  adjustable  0  to  20  dB 


All  active  circuit  elements  are 
solid-state  devices 
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AMPLITUDE 


Fig.  3.5.  BLOCK  DIAGRAM  OF  AMPLITUDE  AND  PHASE- SAMPLING  EQUIPMENT. 

Parameters 


Electronic  Commutator 

2-pole  12-position  switch 
Scans  array  at  100  scans  (over 
12  elements)/sec 
Settling  time:  «  0.3  msec 

100-kHz  Bandpass  Filters 
Bandwidth:  30  kHz 

Phase-Angle  Computer  (Miltron 

Model  351 ) 

Measures  phase  to  1°  or  better 
over  60  dB  signal  variation 
in  either  or  both  channels 

Lowpass  Filters 

Active  2-pole  filter  using  RC 
feedback  around  an  opera¬ 
tional  amplifier 
RC  feedback  set  to  optimize 
rise  time  at  switch  rate 
used 

dc  offset:  <  1  mV 


Linear  Amplitude  Detector 

Linear  to  1%  of  full  scale  over 
dynamic  range  >  AO  dB 


dc  Amplifiers 

dc  offset:  <  1  mV 
Feedback  around  operational 
amplifier  for  good  linearity 


Frequency  Standard  (James  Kniqht 
FS-m  OPT 

100  kHz  output 

Aging  rate  of  1  to  2  parts  in 
1010/day  (1  to  2-1 0-5  Hz/da y 
at  100  kHz) 
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linear  amplitude  detectors  one  of  which  is  used  exclusively  for  ampli¬ 
tude  calibration.  The  detected  amplitude  and  relative  phase  data  samples 
are  low-pass  filtered  and  made  available  for  continuous  monitoring  on  an 
oscilloscope  before  being  processed  for  data  recording. 

A  block  diagram  of  the  data  collection  equipment  is  shown  in  Figure 
3.6.  The  detected  amplitude  and  phase  data  samples  are  multiplexed  and 
converted  to  12-bit  digital  words  under  the  control  of  a  general  purpose 
PDP-8  computer.  The  12-bit  digital  format  allows  for  a  maximum  of  4096 
quantization  levels  in  the  analog-to-digital  (A/D)  conversion  process. 


advance  pulse 


Fig.  3.6.  BLOCK  DIAGRAM  OF  DATA-COLLECTIQN  EQUIPMENT. 


Parameters 


A/D  Converter  and  Multiplexer  Computer  (DEC  PDP-8)  with  Hardware 

(DEC  AFOTT  Arithmetic  Unit 


Multiplexer  switching  speed: 
<  2  ysec 

A-D  word  length:  12  bits 
(3096  levels) 

35  ysec  conversion  time 
+0.25%  max  switching  point 
error  for  +1/2  least  sig¬ 
nificant  bit  of  quantiza¬ 
tion  error 


1.5  psec  memory  cycle  time 
4096  word  fast  memory 
12-bit  word  length 

Magnetic  Tape  Units 

DEC  TU-20  (IBM  Compatible) 
DEC  TU-55  (Program  DECtape) 
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Fig.  3.7.  BLOCK  DIAGRAM  OF  TRANSMITTER 


The  buffer  amplifiers  preceding  the  A/D  converter  are  adjusted  so  that 
the  full  dynamic  range  of  the  converter  is  utilized.  In  converting  the 
phase  measurements,  a  range  of  -  180°  to  +  180°  corresponds  respectively 
to  converter  readouts  of  0  to  12  bits  (0  to  4095  levels);  in  converting 
the  amplitude  measurements,  the  gains  are  set  so  that  no  input  signal 
corresponds  to  0  bits  (0  level)  and  a  fixed  calibration  signal  corre¬ 
sponds  to  11  bits  (2048  level).  After  analog-to-digital  conversion,  the 
amplitude  and  phase  data  samples  are  written  and  stored  on  magnetic  tape. 


3.2.3  Transmitter  Assembly 

A  block  diagram  of  the  transmitting  system  is  shown  in  Figure  3.7. 
The  transmitter  is  essentially  the  same  as  that  used  previously  by  Cox 
(1967),  Cox  and  Waterman  (1971),  Cianos  (1971  and  1978),  Cianos  and 
Waterman  (1973),  and  Waterman  (1973).  The  transportable  transmitting 
antenna  shown  in  Figure  3.8  uses  the  same  2.4-m  (8-ft)  parabolic  reflec¬ 
tor  as  before  but  employs  a  horizontally  polarized  feed.  The  transmit¬ 
ter  has  available  a  maximum  output  power  of  20  watts;  however,  during 
the  actual  tests,  the  nominal  operation  for  the  transmitter  was  between 
10  and  12  watts.  Some  modifications  were  made  in  the  output  power- 
monitoring  assembly  so  that  power  level  adjustments  could  be  performed 
while  the  transmitter  was  in  a  non-radiating  mode.  These  modifications 
included  the  use  of  a  dummy  load  and  additional  coaxial  switches  to 
permit  the  transmitter  signal  to  be  monitored  under  loaded  non-radi¬ 
ating  conditions  during  calibration  and  alignment  periods. 


Fig.  3.8.  TRANSMITTING  ANTENNA  AND  EQUIPMENT  VAN. 


3.3  Experimental  Procedures 

3.3.1  Data  Collection 

Data  was  collected  and  recorded  during  a  continuous  42-hour  experi¬ 
ment  conducted  between  0000  hours  on  21  April  1978  and  1830  hours  on  22 
April  1978.  Ten  minutes  of  transhorizon  data  were  recorded  every  hour 
beginning  approximately  on  the  hour.  The  sampled  amplitude  and  phase 
data  from  each  array  element  were  recorded  in  12-bit  digital  format  on 
magnetic  tape  at  a  40-Hz  rate.  A  five  millisecond  scan  period  was  used 
to  sample  and  record  the  amplitude  and  phase  data  from  all  twelve  ele¬ 
ments.  The  data  samples  were  recorded  In  the  sequence  a  ,$  ;  a  ,<p  ;  .  . 

11  2  2 

.  .  a  . 

12  12 

A  calibration  signal  transmitted  from  a  1 ine-of-sight  source  about 
a  mile  from  the  receiver  along  the  path  was  recorded  for  1.6  seconds  (64 
scans)  before  every  ten-minute  data  recording  period.  Immediately  fol¬ 
lowing  the  recording  of  transhorizon  data,  the  transmitter  was  turned  off 
and  the  calibration  signal  again  turned  on.  The  amplitude  and  phase  cal¬ 
ibration  levels  of  each  array  element  were  then  monitored  with  the  digi¬ 
tal  voltmeter  (Figure  3.5)  and  read  into  the  log  while  the  commutator  was 
manually  stepped  through  all  twelve  element  positions. 

During  actual  data  recording,  the  sampled  amplitudes  and  phases  of 
all  twelve  elements  were  continuously  monitored  on  an  oscilloscope  (Fig¬ 
ure  3.5).  Any  signal  anomalies  that  could  be  attributed  to  aircraft, 
interference,  etc.  were  noted  in  the  log.  Visual  sightings  of  aircraft 
In  the  vicinity  of  the  common  volume  were  logged  at  the  transmitter, 
receiver,  and  PIBAL  launch  sites. 
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3.3.2  System  Calibration 

Prior  to  each  transhorizon  data  recording  period,  the  transmitter 
was  completely  checked  and  aligned;  and  all  array  element  receiving  chan¬ 
nels  were  calibrated.  The  transmitter  frequency  standard  (Figure  3.5) 
was  adjusted  each  day  during  the  quiet  early  morning  hours  and  set  to 
within  2  parts  in  1011  (0.064  Hz  at  3.2  GHz)  of  the  frequency  standard 
used  to  generate  the  reference  and  calibration  signals.  In  the  cali¬ 
bration  and  alignment  of  the  individual  receiving  channels,  a  number  of 
critical  checks  were  made  to  ensure  proper  operation  during  the  data- 
gathering  portions  of  the  experiment.  The  mixer  currents,  the  amplifier 
gains,  and  the  phase-lock  loop  status  of  the  receiving  system  were  rou¬ 
tinely  checked  and  if  necessary,  adjusted  before  each  data  recording 
period. 

During  calibration  periods,  the  transhorizon  signal  was  turned  off 
and  a  3.2001 -GHz  calibration  signal  was  transmitted  from  a  stable  source 
about  1.6  kilometers  from  the  receiver.  This  signal  was  remotely  con¬ 
trolled  from  the  receiving  site  and  was  used  to  properly  set  the  vari¬ 
able  100-kHz  amplifiers  and  phase-shifters  in  each  receiving  channel. 

With  the  calibration  signal  on,  the  commutator  was  manually  advanced 
through  all  twelve  element  positions  and  the  detected  amplitude  and 
adjacent  element  phase-difference  signals  were  monitored  with  the  digi¬ 
tal  voltmeter.  The  variable  phase-shifters  and  amplifiers  were  then 
adjusted  on  an  element  by  element  basis  for  predetermined  calibration 
readings  on  the  digital  voltmeter. 

The  predetermined  phase-difference  and  amplitude  calibration  set- 
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tings  were  obtained  from  the  results  of  a  series  of  1 ine-of-sight  (LOS) 
tests  In  which  the  transmitter  was  placed  at  selected  azimuthal  posi¬ 
tions  about  15  kilometers  from  the  receiving  array.  With  the  calibration 
signal  off  and  the  transmitter  on  (at  reduced  power),  the  individual  am¬ 
plitudes  were  set  uniformly  across  the  array  by  adjusting  the  gains  of 
the  100-kHz  amplifiers.  Adjacent  element  phase-differences  were  set  with 
the  variable  phase-shifters  to  values  calculated  from  the  signal's  known 
angle-of-arrival .  After  the  phase  and  amplitude  adjustments  were  made, 
the  transmitter  was  turned  off  and  the  LOS  calibration  source  turned  on. 
Phase  and  amplitude  readings  were  then  obtained  for  all  twelve  elements 
using  the  digital  voltmeter. 

The  above  procedures  were  repeated  for  three  different  1 ine-of-sight 
transmitter  positions.  Results  from  two  of  the  azimuthal  positions  were 
averaged  and  used  to  obtain  the  final  amplitude  and  phase-difference  cal¬ 
ibration  settings.  The  maximum  variation  in  the  adjacent  element  phase- 
differences  observed  between  the  two  locations  was  14°  and  the  rms  vari¬ 
ation  was  7°.  In  the  amplitude  readings,  the  maximum  variation  was  2  dB 
and  the  rms  variation  was  0.7  dB. 

3.3.3  PIBAL  Wind  Monitoring 

Pilot  balloons  released  into  the  common  volume  from  a  launch  site 
near  Livermore  were  optically  tracked  during  daylight  hours  with  a  single 
theodolite  system  to  obtain  vertical  profiles  of  the  horizontal  wind- 
field.  A  pilot  balloon  (PIBAL)  was  launched  every  hour  at  five  minutes 
past  the  hour  between  0805  and  1805  on  each  of  the  two  data  recording 
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days.  The  balloons  were  released  by  members  of  the  San  Jose  State  Uni- 
verstty  Meteorology  Department  from  a  position  beneath  the  common  volume 
about  0.9°  SE  of  the  great-circle  bearing  and  165  meters  (540  ft)  above 
mean  sea  level.  The  PIBAL  ascent  rate  was  approximately  200  meters  (650 
ft)  per  minute.  Azimuth  and  elevation  position  readings  were  taken  every 
30  seconds  and  each  PIBAL  was  tracked  to  an  altitude  of  approximately 
2.4  km  (7870  ft)  above  sea  level.  An  algorithm  used  by  the  meteorology 
group  reduced  the  PIBAL  azimuth  and  elevation  information  to  horizontal 
wind  velocity  profiles. 
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Chapter  4 


ANALYTICAL  METHODS 

Two  basic  signal  processing  steps  are  used  to  obtain  the  angle-of- 
arrival  Doppler  spectra  from  measured  amplitude  and  phase  data  samples. 
First,  a  digital  Fourier  transform  is  performed  in  the  time  domain  on  the 
sampled  output  signals  from  the  individual  array  elements.  The  ampli¬ 
tudes  and  phases  of  the  resultant  spectra  are  preserved  and  used  to 
characterize  the  array  element  output  signals  in  the  frequency  domain. 

A  spatial  transform  is  then  performed  in  the  array  plane  on  the  element 
spectra  to  obtain  an  angular  Doppler  spectrum  for  the  array.  In  this 
process,  the  individual  element  spectra  are  combined  using  conventional 
additive  array  techniques  to  maximize  the  output  spectra  for  a  particu¬ 
lar  angle-of-arrival .  An  average  Doppler  frequency  shift  (spectrum  cen¬ 
troid)  is  computed  as  a  function  of  the  angular  "pointing"  direction  of 
the  array  with  respect  to  the  great-circle  (zero-Doppler)  plane.  The 
apparent  source  direction  together  with  the  computed  average  Doppler 
shift  is  used  to  infer  an  average  transverse  wind  velocity. 

4.1  Single  Element  Spectra 

The  measured  element  amplitude,  an(k)  ,  and  relative  phase, 

4>n Ck)  >  data  samples  are  used  to  construct  the  single  element  Doppler 
spectra.  Each  spectral  component  is  evaluated  from  a  set  of  K  con¬ 
tiguous  data  samples.  The  in-phase  and  out-of-phase  quadrature  com¬ 
ponents  of  the  element  spectra  are  computed  in  accordance  with  the 
following 


39 


V" >  r  2irf .  (k-1 )  "1 

SXK^n,f1^  =  X  ,an(k)cos  L  40  ^n(k)J 


r  2irf . 

SYK(n-V  =  2^a"(k)Sin  L~ 


2^  (k-1) 
40 


■  -  ♦„(«] 

(4.1) 

--*n'k>] 

(4.2) 

2,  ....  12 

The  amplitudes,  AJf. )  ,  and  phases,  ¥  „(f. )  ,  of  the  spectral  com- 

n  ^  1  n  k  i 

ponents  are  computed  from  the  quadrature  phasor  components 


AnK^fi ^  =  '  SXK^n,fi ^  +  SYK^n,fi ^ 

^nK^f1  ^  =  tan  [SYK^n,fi^/SXK^n,fi?l 


(4.3) 

(4.4) 


These  parameters  are  obtained  for  100  Doppler  frequencies  in  a  preselect¬ 
ed  spectral  observation  window.  The  spectral  window,  which  is  less  than 
the  receiver  bandwidth  (23  Hz),  is  centered  about  zero  Doppler  to  permit 
simultaneous  observation  of  the  upwind  (positive)  and  downwind  (negative) 
Doppler  shifts. 

4.2  Array  Angular  Spectra 

The  array  angular  spectra  are  derived  from  the  amplitudes,  A„„(f. ) 
and  phases,  ¥  ^(tj)  ,  of  the  individual  element  spectra.  A  progressive 

linear  phase-shift  across  the  array  is  applied  to  the  spectral  components 
of  the  individual  element  spectra.  Each  component  is  phase-shifted  in 
accordance  with  conventional  additive  array  techniques.  The  quadrature 
phasor  components  of  the  resultant  angular  spectra  are  then  computed  for 
selected  beam-pointing  angles,  a  ,  between  +  1.94°  as  follows 


(4.5) 


SXK(VV 

SYK{ar’V 
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where:  y  =  phase  shift  constant 

=  92.71  phase  degrees/beam  angle  degree 


A  positive  (negative)  angular  direction  corresponds  to  a  SE  (NW)  angle- 
of-arrival  with  respect  to  the  great-circle  bearing.  The  Doppler  or  po¬ 
wer  spectra,  P(ar,fi )  ,  are  computed  for  a  particular  angular  direction 

by  summing  the  squares  of  the  quadrature  phasor  components,  i.e. 

P(ar,f.)  =  SjK(ctr,f.)  +  Sj^.f.)  (4.7) 


These  angular  spectra  are  essentially  unweighted  since  each  spectral  com¬ 
ponent  is  evaluated  from  the  same  set  of  element  amplitudes  irrespective 
of  its  source  direction.  The  Doppler  spectra,  however,  can  be  weighted 
to  account  for  the  angl e-of-arrival  variations  arising  from  the  scatter¬ 
ing  mechanism  and  the  antenna  patterns.  We  define  the  direction  to  the 
source  in  terms  of  an  angle,  a..  ,  viewed  from  the  receiver.  The 

weighted  Doppler  spectra,  P  (a  ,f. )  ,  are  then  determined  by  applying 

W  I  • 

a  normalized  angular  scattering  function,  W^(a^)  ,  and  a  combined  re¬ 
ceiver-transmitter  antenna  gain  function,  Gc(a.-ar)  ,  to  the  unweighted 
spectra,  i.e. 

Pw(ar,fi)  =  P(ar,fi)  Wg(a.)  G^ot.-oJ  (4.8) 
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Each  Doppler  frequency  is  converted  to  an  apparent  angular  source  direc¬ 
tion  through  use  of  the  Doppler  relation  In  (2.15),  which  for  a  uniform 
transverse  wind  velocity*,  vQ  ,  yields 


(4.9) 


On  the  basis  of  the  computed  angular  source  direction,  each  spectral  com¬ 
ponent  Is  weighted  in  accordance  with  the  classical  isotropic  turbulence 
model,  a  sin(x)/x  receiving  array  pattern,  and  a  Gaussian  transmitting 
pattern.  The  normalized  scattering  function  Is  evaluated  for  a  given  ele¬ 
vation  angle,  3  ,  using  the  following  expression  with  the  numerator  ob¬ 

tained  from  equation  2.20 


M3(o1)  = 


1.51  +  |  V  a2.  +  32 


_-m 


0.84 


(4.10) 


where:  m  =  11/3 
8  >  0.25° 

The  combined  antenna  gain  for  the  assumed  patterns  is  evaluated  for  each 
source  direction  from  the  following 


G  (a. -a  ) 
cv  l  r 


sin[a(ai-ar)] 


air(a.-a  )/l80 
i  r 


exp[-  ba?] 


(4.11) 


where:  a  =  549.900 
b  =  0.03429 


The  sign  convention  for  the  wind  velocity  is  negative  for  a  SE  wind  and 
positive  for  a  NW  wind. 
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The  beam  constants  a  and  b  are  determined  for  half-power  beamwidths 
of  0.29°  and  3°,  respectively.  The  transmitting  beam  constant,  b  , 
includes  the  path  geometry  scaling  factor  which  permits  angular  measure¬ 
ments  to  be  referenced  to  the  receiver. 

4.3  Average  Doppler  Shift 

The  average  Doppler  shift  is  obtained  by  computing  the  first  moment 
(spectrum  centroid)  of  the  angular  Doppler  spectra.  For  the  unweighted 
spectra,  the  average  is  given  by 
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100 

£p'vfi' 
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Similarly,  for  the  weighted  spectra 


(4.12) 
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4.4  Average  Wind  Velocity 

The  average  transverse  wind  velocity  is  determined  using  the  results 
of  averaging  the  unweighted  and  weighted  Doppler  spectra.  An  unweighted 
average  velocity,  v  ,  is  obtained  directly  from  the  measured  unweighted 
average  Doppler  shift,  f  ,  by  using  equation  2.15,  i.e. 
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A  weighted  average  wind  velocity,  v  ,  is  computed  to  Include  correc- 
tlons  for  the  effects  of  the  scattering  mechanism  and  the  antenna  pat¬ 
terns.  This  result  is  obtained  from  a  linear  Interpolation  of  the  mea¬ 
sured  weighted  average  Doppler  shift,  f  ,  (for  uniform  wind  v  )  and 

w  o 

the  measured  unweighted  average  Doppler  shift 


(4.15) 


In  these  expressions,  the  wind  speed  is  measured  in  meters  per  second 
and  all  angles  are  measured  in  degrees. 


Chapter  5 


DISCUSSION  OF  DATA 


5.1  PIBAL  Wind  Data 

Local  wind  conditions  were  influenced  during  the  tests  by  a  low 
pressure  weather  system  centered  off  the  coast  of  the  state  of  Washing¬ 
ton.  About  six  hours  before  the  onset  (0100  April  21)  of  the  experiment, 
a  low  pressure  frontal  system  moved  through  the  area  with  slight  preci¬ 
pitation  and  some  gusty  winds.  During  the  two  11 -hour  periods  between 
0800  and  1800  when  pilot  balloons  (PIBALS)  were  launched,  the  measured 
cross-path  winds  were  generally  from  the  southeast.  Steady  southeast¬ 
erly  winds  prevailed  at  nearly  all  altitudes  on  the  first  day  of  the  ex¬ 
periment.  The  winds  were  moderate  and  progressively  increased  during 
the  day  reaching  maximum  velocities  in  mid-afternoon.  In  contrast, 
light  variable  winds  persisted  throughout  the  second  day  when  wind  speeds 
rarely  exceeded  two  meters  per  second.  Typical  transverse  winds  on  this 
day  were  less  than  one  meter  per  second  and  on  occasions  wind  direction 
changes  were  observed  between  SE  and  NW. 

PIBAL  readings  were  taken  every  30  seconds  at  approximately  100-m 
height  intervals.  PIBAL  data  above  the  radio  horizon  were  averaged  over 
400-m  vertical  segments  for  comparison  with  inferred  results  from  radio 
Doppler  measurements.  Successive  averages  were  taken  throughout  the 
vertical  extent  (1600  meters)  of  the  common  volume  region.  Average 
transverse  wind  velocities  in  the  four  vertical  segments  of  the  common 
volume  are  tabulated  in  Tables  5.1  and  5.2  for  the  PIBAL  wind  soundings 
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Table  5.1 

PIBAl  HIND-MEASUREMENT  QATA  ON  21  APRIL  1978* 


PIBAL  Average  Transverse  Wind  Velocity  (m/s) 


Release 

1st  400-m 

2nd  400-m 

3rd  400-m 

4th  400-m 

0805 

1.93  SE 

3.58  SE 

5.69  SE 

7.12  SE 

0905 

1.06  SE 

3.37  SE 

5.61  SE 

6.85  SE 

1005 

1.16  SE 

2.73  SE 

5.19  SE 

5.86  SE 

1105 

1 .84  SE 

1  .96  SE 

3.96  SE 

5.12  SE 

1205 

4.10  SE 

4.49  SE 

1 .68  SE 

4.19  SE 

1305 

3.95  SE 

4.18  SE 

4.55  SE 

4.02  SE 

1405 

3.21  SE 

3.96  SE 

4.81  SE 

4.01  SE 

1505 

1 .98  SE 

1.08  SE 

3.51  SE 

4.51  SE 

1605 

3.52  SE 

4.32  SE 

4.74  SE 

4.65  SE 

1705 

4.96  SE 

4.66  SE 

2.51  SE 

3.55  SE 

1805 

3.15  SE 

2.40  SE 

2.92  SE 

3.69  SE 

Table  5.2 

PIBAL  HIND-MEASUREMENT  DATA  ON  22  APRIL  1978* 


PIBAL  Average  Transverse  Wind  Velocity  (m/s) 


Release 

1st  400-m 

2nd  400-m 

3rd  400-m 

4th  400-m 

0805 

0.10  NW 

0.92  NW 

0.48  SE 

2.82  SE 

0905 

0.12  NU 

1.49  NW 

0.14  SE 

1.76  SE 

1005 

0.34  SE 

0.09  SE 

0.53  NW 

0.98  SE 

1105 

0.54  NW 

0.14  NW 

1.00  NW 

0.08  SE 

1205 

0.43  NU 

0.78  SE 

0.99  NW 

1.13  NW 

1305 

0.09  SE 

0.29  SE 

2.17  NW 

0.40  NW 

1405 

0.53  SE 

1.88  SE 

0.56  NW 

2.25  NW 

1505 

1.45  SE 

0.73  SE 

1.07  NW 

0.09  NW 

1605 

1.01  SE 

0.72  SE 

0.98  NW 

0.52  NW 

1705 

0.45  SE 

0.22  SE 

0.70  NW 

2.72  NW 

1805 

0.89  NW 

0.18  NW 

2.17  NW 

1.02  NW 

*Average  transverse  wind  in  400-m  vertical  segments  above  radio  horizon. 
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taken  respectively  on  April  21  and  22.  Examples  of  vertical  profiles  of 
the  cross-path  wind  component  observed  on  the  two  days  are  shown  in  Fig¬ 
ure  5.1  [(a)  thru  ( f ) ] .  Data  from  the  0805  launch  on  April  21  (Figure 
5.1a)  reflected  a  fairly  strong  wind  shear  between  the  1300  and  1600- 
meter  elevations.  Within  this  300-meter  vertical  range,  the  wind  changed 
from  0.8  m/s  NW  to  12  m/s  SE.  This  shear  dissipated  by  the  next  wind 
sounding  one  hour  later.  No  other  anomalous  wind  conditions  were  ob¬ 
served  for  the  duration  of  the  experiment. 

5.2  Wind-Produced  Doppler 

The  Doppler  spectra  presented  in  this  section  are  computed  from  600 
contiguous  sets  of  data  samples  representing  15  seconds  of  real-time  data. 
Each  data  set  is  comprised  of  twelve  amplitudes  and  twelve  relative 
phases.  Computations  are  made  from  the  same  data  sets  for  a  sequence  of 
array  pointing  angles  about  the  great-circle  bearing.  In  order  to  simul¬ 
taneously  measure  upwind  and  downwind  Doppler  shifts,  spectral  observa¬ 
tion  windows  are  centered  about  the  zero-Doppler  frequency.  Observation 
windows  used  in  the  spectrum  analysis  are  narrower  than  the  system  noise 
bandwidth  (23  Hz)  but  are  broad  enough  to  encompass  the  ambient  wind- 
produced  Doppler  shifts.  For  the  wind  conditions  experienced  during  the 
tests,  an  8  to  10-Hz  wide  observation  window  was  found  to  be  sufficient 
to  completely  envelop  the  wind-produced  Doppler  spectrum.  Outside  the 
selected  observation  windows,  the  Doppler  spectra  have  relatively  low 
power  spectral  densities.  Each  computed  spectrum  contains  100  Doppler 
frequencies  spaced  from  about  0.08  to  0.10  Hz  apart  depending  upon  the 
width  of  the  spectral  observation  window.  The  spectral  resolution 


Fig.  5.1.  CONTINUED 
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attainable  in  the  coherent  integration  process  of  the  600  data  samples 
is  0.07  Hz.  However,  the  realizable  resolution  is  more  on  the  order  of 
0.1  Hz  due  to  the  short-term  instabilities  (1.5  parts  in  lO^1)  of  the 
quartz  oscillators. 

5.2.1  Unweighted  Angular  Spectra 

Examples  of  some  unweighted  Doppler  spectra  are  shown  in  Figure  5.2 
[(a)  thru  (1)].  These  spectra  are  taken  for  symmetrical  (NW/SE )  beam 
pointing  directions  about  the  great-circle  bearing  at  0°  azimuth.  The 
spectra  are  plotted  over  a  40-dB  range  in  0.25°  increments  out  to  one 
degree  either  side  of  the  great-circle  bearing.  Two  15-second  samples 
taken  one  minute  apart  are  presented  for  each  hourly  example  in  order  to 
show  temporal  variations  in  the  wind-produced  Doppler  shifts. 

Initial  observations  of  the  high-resolution  spectra  show  the  pres¬ 
ence  of  a  coherent  stationary  source  with  a  very  narrow  spectral  width 
at  zero-Doppler  frequency.  This  spectral  component  is  most  likely  a  re¬ 
sult  of  knife-edge  diffraction  over  one  or  more  of  the  mountain  ridges 
between  the  transmitter  and  receiver.  It  is  seen  that  the  diffracted 
component  increases  in  intensity  as  the  array  scan  approaches  the  great- 
circle.  Occasionally,  when  the  array  beam  is  pointed  on-axis,  the  level 
of  the  diffracted  component  exceeds  the  40-dB  plotting  range  and  a  lin¬ 
ear  foldover  (overshoot)  occurs.  Aside  from  its  intensity  level,  the 
diffraction  component  is  totally  invariant  to  the  azimuthal  pointing 
direction  of  the  array. 

In  contrast,  the  wind-produced  Doppler  shifts  show  systematic  varia¬ 
tions  with  the  array  pointing  direction.  A  qualitative  analysis  of  the 
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Fig.  5.2.  CONTINUED. 
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(d) 

Mg.  5.2.  CONTINUED. 
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Fig.  5.2.  CONTINUED. 
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Fig.  5.2.  CONTINUED. 
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Fig.  5.2.  CONTINUED. 
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wind-produced  spectra  shows  progressive  increases  in  the  average  Doppler 
shift  as  the  array  is  pointed  off-axis.  For  the  examples  shown  in  Figure 
5.2,  SE  winds  are  apparent  at  all  angles.  Doppler  shifts  are  generally 
positive  for  array  pointing  angles  in  the  SE  (upwind)  direction  and  nega¬ 
tive  for  angles  in  the  NW  (downwind)  direction.  At  most  off-axis  angles, 
the  upwind  and  downwind  Doppler  spectra  are  nearly  antisymmetric  about 
the  zero-Doppler  frequency  which  suggests  some  degree  of  spatial  uni¬ 
formity  in  the  transverse  winds.  On-axis  the  Doppler  spectra  are  gen¬ 
erally  symmetrical  about  zero  Doppler.  Any  asymmetries  appearing  in  the 
spectra  when  the  array  beam  is  aligned  along  the  great-circle  can  either 
be  attributed  to  a  vertical  wind  component  or  to  spatial  non-uniformities 
in  the  transverse  wind.  Birkemeier  et  al  (1977)  point  out  that  such 
spectral  asymmetries  also  result  from  the  misalignment  of  the  trans¬ 
mitting  and  receiving  antennas.  Misalignment  errors,  however,  are  elimi¬ 
nated  when  using  a  phased-array. 

Although  no  quantitative  analysis  is  presented,  other  observations 
from  Figure  5.2  reveal  definite  spectral  broadening  during  afternoon 
periods  when  the  atmosphere  is  fairly  well  mixed.  This  is  a  logical 
consequence  of  diffuse  (turbulent)  scattering  processes  which  produce 
wide  variations  in  the  signal's  angle-of-arrival .  On  the  other  hand, 
in  the  morning  hours,  the  atmosphere  is  relatively  stable  and  is  more 
likely  to  be  partially  stratified.  In  such  situations,  a  quasi-coherent 
signal  is  received  that  arrives  via  a  number  of  specular  reflection 
points  distributed  over  a  small  angular  region  near  the  great-circle 
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bearing.  As  a  result  of  the  narrow  range  of  angles-of-arrival ,  the  sig¬ 
nal  spectra  are  not  nearly  as  broad  as  the  spectra  produced  by  signals 
propagating  through  a  turbulent  atmosphere.  The  spectra  are  also  seen 
to  be  somewhat  broader  on  the  first  day  [Figure  5.2  (a)  thru  (f)]  of  the 
experiment  than  on  the  second  day  [Figure  5.2  (g)  thru  (1)].  This  vari¬ 
ance  appears  directly  attributable  to  the  contrasting  wind  conditions  on 
the  two  days. 

The  short-term  variability  of  the  wind  can  be  observed  from  temp¬ 
oral  changes  in  the  Doppler  spectra.  Of  particular  importance  in  this 
regard  are  the  temporal  changes  in  the  average  Doppler  shift  as  a  func¬ 
tion  of  the  array  pointing  angle.  Figure  5.3  [(a)  thru  (f)]  shows  the 
measured  average  Doppler  shifts  for  the  spectra  given  in  Figure  5.2. 
Average  Doppler  shifts  are  computed  in  accordance  with  equation  4.12  for 
21  equally-spaced  (0.15°)  pointing  angles  between  1.5°  NW  and  1.5°  SE. 
The  slopes  of  these  curves  are  directly  proportional  to  the  prevailing 
wind  speeds.  Averages  taken  one  minute  apart  are  shown  on  the  same 
plot  for  comparative  analysis  of  short-term  variations.  The  15-second 
samples  are  generally  in  close  agreement  from  one  minute  to  the  next. 
This  suggests  reasonably  good  correlation  of  wind  conditions  within  one 
minute  intervals.  As  evidenced  by  the  near  linearity  of  the  Doppler 
shift  curves,  there  is  some  apparent  uniformity  in  the  wind  structure 
over  the  spatial  regions  probed.  Occasionally,  however,  there  are 
noticeable  departures  from  this  linear  trend  that  may  be  attributable 
to  variable  wind  conditions.  For  example,  in  Figure  5.3a  at  0803  on  21 
April  a  definite  change  in  the  polarity  of  the  average  Doppler  shift  is 
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g.  5.3.  MEASURED  UNWEIGHTED  AVERAGE  DOPPLER  SHIFT  VS 
ARRAY  POINTING  ANGLE. 
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Fig.  5.3.  CONTINUED. 
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observed  when  the  array  is  pointed  in  the  vicinity  of  0.3°  NW.  An  off- 
axis  polarity  change  such  as  this  is  reflective  of  a  corresponding  wind 
direction  change.  Supporting  this  are  the  wind-shear  conditions  re¬ 
flected  in  the  PIBAL  measurement  data  (Figure  5.1a)  at  the  same  hour. 
There  are  nonlinearities  in  the  average  Doppler  shift  curves  that  are 
not  necessarily  associated  with  the  wind.  These  nonlinearities  occur 
near  the  extremities  of  the  array  sector  scan  and  are  most  likely  due 
to  the  combined  off-axis  effects  of  atmospheric  scattering  and  secondary 
grating  lobes.  Effects  of  the  scattering  mechanism  can  be  taken  into 
account  by  weighting  the  Doppler  spectra  in  accordance  with  an  appropri¬ 
ate  atmospheric  scattering  model.  The  data  results  from  this  effort  are 
discussed  in  the  following  section. 


5.2.2  Weighted  Angular  Spectra 


The  Doppler  spectra  presented  in  the  previous  section  are  essen¬ 
tially  unweighted  in  that  the  spectral  components  are  derived  from  data 
samples  that  are  uniformly  weighted  over  the  angular  range  scanned  by 
the  array.  Here  we  take  into  account  the  apparent  source  direction  of 
each  spectral  component  and  the  angular  scattering  properties  of  the 
atmosphere.  Weighting  factors  based  on  angle-of-arrival  scattering 
variations  are  then  applied  to  the  individual  spectral  components  in  ac¬ 
cordance  with  procedures  outlined  in  section  4.2. 

Figure  5.4  [(a)  thru  (1)3  shows  the  results  of  applying  a  normal¬ 
ized  Tatarski  scattering  model  to  the  unweighted  spectral  data  in  Fig¬ 
ure  5.2.  These  spectra  are  obtained  for  an  assumed  uniform  wind  and  a 
turbulent  scattering  region  situated  at  an  elevation  angle  of  0.25°. 
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g.  5.4.  WEIGHTED  WIND-PRODUCED  DOPPLER  SPECTRA.  Angle-of-arrlval  spec¬ 
tra  are  weighted  In  accordance  with  (1)  a  normalized  Tatarskl  scattering 
model,  (2)  a  s1n(x)/x  receiving  array  pattern,  and  (3)  a  Gaussian  trans¬ 
mitting  pattern  under  the  assumption  of  a  uniform  wind  and  a  source  ele¬ 
vation  angle,  B  ■  0.25°. 


RELATIVE  POWER  (dB) 


RELATIVE  POWER 


(f) 

Fig.  5.4.  CONTINUED. 
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ABSTRACT 

Conventional  methods  for  monitoring  winds  aloft  have  inherent  spa¬ 
tial  and/or  temporal  limitations  that  severely  restrict  their  useful¬ 
ness.  An  alternate  solution  is  the  use  of  transhorizon  bistatic  radio 
techniques.  By  using  a  phase-coherent  system  with  some  form  of  beam¬ 
swinging  capability,  it  is  possible  to  measure  the  wind-produced  Doppler 
shifts  associated  with  moving  off-path  scatterers. 

In  this  experiment  a  rapid  data-gathering  system  is  used  to  re¬ 
motely  probe  the  wind-field  velocity  structure  in  the  common  volume  re¬ 
gion  of  a  164-km  transhorizon  propagation  path.  The  phase-stable  sys¬ 
tem  operates  at  3.2  GHz  and  employs  a  12-element  horizontal  receiving 
antenna  array.  The  amplitudes  and  relative  phases  of  the  output  sig¬ 
nals  from  the  individual  elements  are  sequentially  sampled  and  recorded 
4  at  a  40-Hz  rate.  Off-line  processing  of  the  amplitude  and  phase  data  /.  i 
samples  is  used  to  *scan°  in  azimuth  a  0.3°  wide  fan  beam  through  a 
3.88®  sector  and  to  obtain  virtual  simultaneous  Doppler  spectra  measure¬ 
ments  at  a  continuum  of  angles  within  the  sector  scanned.  An  algorithm 
developed  to  identify  the  apparent  source  direction  of  each  spectral 
component  permits  an  appropriate  weighting  to  be  applied  to  the  spectra 
on  the  basis  of  angl e-of-arrival .  The  average  Doppler  shift  (spectrum 
centroid)  as  a  function  of  the  array  pointing  angle  is  interpreted  in 
terms  of  the  average  transverse  wind  speed  and  direction.  Doppler  re¬ 
lations  applicable  to  previous  experiments  are  modified  fo\account  for 
the  influence  of  differences  in  the  longitudinal  path  geometry.  An 
analytical  expression  employing  a  path  asymmetry  factor  is  developed 
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for  the  wind-produced  Doppler  from  a  single  scatterer  model.  Wind 
velocities  inferred  from  the  average  Doppler  shift  data  are  compared 
with  in  situ  wind  observation  data  obtained  by  optically  tracking  pi¬ 
lot  balloons  released  into  the  common  volume. 

.  Measurement  results  show  the  average  Doppler  shifts  to  have  an 
almost  linear  dependence  on  angl e-of-arrival .  Wind-produced  Doppler 
shifts  resulting  from  average  transverse  velocities  up  to  about  5  m/s 
have  been  measured  and  found  consistent  both  in  magnitude  and  direc¬ 
tion  with  the  prevailing  winds  observed  in  the  lower  strata  of  the 
common  volume  region.  On  numerous  occasions  Doppler  measurement  re¬ 
sults  reflected  much  higher  average  velocities  ('40  m/s)  than  could 
be  attributed  to  normal  low-altitude  winds.  Time-spatial  analyses 
of  these  anomalous  spectra  are  indicative  of  aircraft-induced  atmos¬ 
pheric  disturbances.  The  spectral  signatures  of  these  disturbances 
may  provide  additional  insight  into  the  generation  and  decay  of  air¬ 
craft  trailing  vortices. 
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Chapter  1 


INTRODUCTION 

This  dissertation  describes  the  results  of  a  microwave  transhorizon 
experiment  in  which  a  rapid  data-gather ing  antenna  array  is  used  as  a 
remote  sensing  probe  of  low-altitude  tropospheric  winds. 

1 .1  Background 

The  initial  transhorizon  propagation  of  an  ultrahigh  frequency  sig¬ 
nal  was  reported  by  Marconi  in  1932  after  he  successfully  transmitted  a 
500-MHz  signal  over  a  distance  of  168  miles  (270  km).  The  underlying 
physical  principles  involved  in  transhorizon  propagation  subsequently 
became  a  subject  of  discussion  and  some  controversy  among  researchers. 
The  smooth-earth  diffraction  theory  (Watson,  1918  and  1919;  Van  der  Pol 
and  Bremmer,  1937,  1938,  and  1939),  which  suggests  an  exponential  de¬ 
crease  in  the  received  power  with  distance,  was  the  generally  accepted 
explanation  for  transhorizon  propagation  in  the  early  1940s.  With  the 
advent  of  high-power  transmitters  and  sensitive  receivers  during  the 
war,  signal  reception  far  beyond  the  horizon  became  more  frequent  (Bul- 
lington,  1955).  The  received  signal  levels,  however,  were  consistently 
found  to  be  far  above  the  levels  predicted  by  the  diffraction  theories. 
Concerted  efforts  were  made  by  investigators  in  the  late  1940s  to  find 
a  plausible  explanation  for  the  transhorizon  signal  observations.  Pe- 
keris  (1947)  suggested  the  signals,  based  on  their  appearance,  were  a 
result  of  some  scattering  process  from  refractive-index  variations  in 
the  atmosphere.  As  a  consequence,  several  theoretical  propagation 


1 


models  evolved  that  took  into  account  the  prevailing  atmospheric  struc¬ 
ture.  The  preponderance  of  experimental  data  taken  over  the  years  sup¬ 
ports  the  validity  of  two  basic  theoretical  models,  one  derived  from  a 
turbulence  scattering  theory  (Booker  and  Gordon,  1950;  Staras,  1952) 
and  the  other  derived  from  a  layer  reflection  theory  (Bauer,  1956;  Friis 
et  al.,  1957): 

(1)  In  the  turbulence  theory  the  propagation  phenomena 
is  characterized  as  a  scattering  process  resulting 
from  turbulent  conditions  in  the  atmosphere  that 
cause  random  irregularities  in  the  refractive  pro¬ 
perties  of  the  medium. 

(2)  In  the  layer  theory  the  propagation  phenomena  is 
characterized  as  a  reflection  process  resulting 
from  stratified  layers  in  the  atmosphere  that  pro¬ 
duce  sharp  discontinuities  in  the  refractive  index 
gradient. 

Much  of  the  experimental  and  theoretical  work  conducted  since  the 
introduction  of  the  two  fundamental  theories  has  dealt  with  the  study 
of  the  physical  mechanisms  supporting  propagation  beyond  the  horizon. 

The  early  experimental  research,  however,  was  influenced  to  a  large  ex¬ 
tent  oy  the  demands  for  reliable  long-distance  radio  communications 
(Bullington,  1950;  Mellen  et  al.,  1955).  Consequently,  one  of  the 
principal  concerns  of  researchers  was  the  signal  fading  problems  often 
experienced  on  transhorizon  paths.  The  cumulative  results  of  the  nu¬ 
merous  experiments  conducted  up  to  about  1960  provided  sufficient  data 
to  establish  a  statistical  basis  for  predicting  (Baris  et  al . ,  1962) 
the  average  expected  level  of  performance  of  transhorizon  radio  trans¬ 
mission  systems. 

By  1960  the  tropospheric  propagation  characteristics  beyond  the 
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horizon  were  fairly  well  understood.  Investigators  subsequently  began 
to  look  into  ways  of  distinguishing  between  the  propagation  mechanisms 
associated  with  the  turbulence  and  layer  model  theories.  At  Stanford 
researchers  (Waterman  et  al.,  1957;  Waterman,  1958;  Lee,  1961;  Stroh- 
behn,  1963;  Waterman  and  Strohbehn,  1963;  Strohbehn  and  Waterman,  1966; 
Cox,  1967;  Cox  and  Waterman,  1969  and  1971;  Cianos,  1971;  Cianos  and 
Waterman,  1973;  Waterman,  1973;  Cianos,  1978)  have  used  various  experi¬ 
mental  techniques  to  differentiate  between  the  two  prevailing  atmos¬ 
pheric  structures  suggested  by  the  theoretical  models.  These  empirical 
methods  provided  in  effect  a  means  of  remotely  probing  the  atmosphere 
with  radio  waves  —  a  technique  in  which  analyses  of  specific  parameters 
of  the  received  signal  are  used  to  infer  certain  characteristics  about 
the  intervening  medium. 

Several  important  meteorological  parameters  can  be  measured  re¬ 
motely  using  transhorizon  radio  techniques.  Among  the  many  measureable 
parameters,  the  wind  velocity  structure  in  the  lower  atmosphere  has  re¬ 
ceived  considerable  attention  by  investigators  (Atlas  et  al.,  1968; 
Birkemeier  et  al . ,  1968;  Atlas  et  al . ,  1969;  Atlas,  1969;  Birkemeier  et 
al.,  1969;  Lammers  and  Olsen,  1973;  Chadwick  et  al . ,  1975)  over  the 
past  decade.  Early  investigations  of  the  signal  fading  phenomena  on 
transhorizon  paths  have  demonstrated  marked  correlations  between  signal 
fade  rates  and  local  transverse  wind  conditions  (Laaspere,  1958;  Craw¬ 
ford  et  al . ,  1959;  Doherty,  1959;  Gjessing,  1964).  Refractive  struc¬ 
tures  within  the  common  volume  drift  with  the  mean  wind  velocity  which 
results  in  stochastic  phase  variations  in  the  received  signal.  These 
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phase  variations  are  reflected  in  the  Doppler  spectrum  and  are  depen¬ 
dent  upon  the  spatial  distribution,  size,  and  relative  velocities  of 
the  moving  atmospheric  structures.  The  analytical  and  experimental 
work  performed  in  this  area  have  provided  the  motivation  for  the  ex¬ 
periment  described  in  this  dissertation  which  seeks  to  further  our  un¬ 
derstanding  of  the  spatial  and  temporal  variations  of  atmospheric  mo¬ 
tion  produced  by  the  wind. 

1 .2  Significance  of  Experiment 

Knowledge  of  wind  information  in  localized  regions  of  the  lower 
atmosphere  is  of  vital  importance  not  only  to  meteorologists  but  to 
researchers  in  such  diverse  areas  as  air-traffic  safety,  environmental 
pollution,  energy,  and  even  agriculture.  The  conventional  rawinsonde 
technique  of  monitoring  winds  aloft  has  inherent  spatial  and  temporal 
limitations  that  severely  restrict  its  usefulness.  Situations  requir¬ 
ing  near  real-time  wind  information  in  selected  regions  of  the  atmo¬ 
sphere  are  not  adequately  served  by  the  widely  spaced  and  infrequent 
(every  12  hours)  launches  of  rawinsondes.  A  reduction  in  either  the 
time  interval  or  the  spacing  between  launches  becomes  economically  pro¬ 
hibitive  in  terms  of  the  increased  manpower  and  instrumentation  re¬ 
quirements. 

An  alternate  means  of  wind  sounding  in  the  lower  atmosphere  has 
recently  been  proposed  by  Chadwick  et  al.  (1975).  The  remote  probing 
technique  proposed  involves  the  use  of  a  network  of  bistatic  radio 
transmission  systems.  The  proposal  stems  from  the  results  of  experi¬ 
mental  work  by  Birkemeier  et  al.  (1968  and  1969)  and  Laimers  and  Olsen 
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(1973).  Analyses  by  Atlas  (1969)  and  Atlas  et  al.  (1965,  1968,  and 
1969)  have  provided  considerable  insight  into  the  elementary  theory  and 
have  enabled  the  experimental  observations  to  be  interpreted  in  meteoro¬ 
logically  significant  terms.  Observed  systematic  increases  in  the  aver¬ 
age  Doppler  shift  with  offset  angle  from  the  great-circle  bearing  have 
been  attributed  to  the  wind-produced  drift  of  refractive  structures  in 
the  atmosphere.  By  incorporating  narrow-beam  antennas  with  some  form  of 
beam-swinging  capability,  it  is  possible  to  remotely  monitor  the  trans¬ 
verse  (cross-path)  component  of  the  wind  over  a  wide  spatial  extent  using 
a  bistatic  phase-coherent  transhorizon  system.  Aligning  two  such  systems 
orthogonal  to  one  another  permits  a  complete  characterization  of  the  hori 
zontal  wind-field. 

Previous  experiments  have  employed  synchronously  offset  steerable 
antennas  to  probe  selected  regions  in  the  atmosphere  and  have  been  con¬ 
ducted  on  longitudinally  symmetrical  propagation  paths.  These  are  trans¬ 
horizon  paths  where  the  common  volume  region  is  symmetrical  about  the  mid 
path  plane  between  the  receiver  and  transmitter  terminals.  Over  highly 
asymmetrical  paths  the  common  volume  region  is  skewed  significantly  to 
one  side  of  the  mid-path  plane.  The  resultant  Doppler  relations  must 
then  reflect  the  influence  of  the  longitudinal  position  of  the  dominant 
scattering  region.  Analytical  results  applicable  to  both  longitudinally 
symmetrical  and  asymmetrical  paths  are  presented  in  this  dissertation 
along  with  experimental  measurement  data.  Unique  features  of  the  experi¬ 
ment  pertinent  to  the  remote  probing  of  tropospheric  winds  are 

(1)  Angl e-of-arrival  measurements  of  wind-produced  Doppler 
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are  derived  from  processing  array-sampled  amplitude 
and  relative  phase  data. 

(2)  Rapid  sampling  features  of  the  data-gathering  array 
permit  the  simultaneous  measurement  of  the  upwind 
and  downwind  Doppler  spectra. 

(3)  Measurements  are  made  on  a  longitudinally  asymmetrical 
propagation  path  where  the  principal  scattering  region 
is  centrally  located  between  the  receiver  and  the  mid¬ 
path  plane. 

1 .3  Description  of  Experiment 

To  investigate  the  spatial  and  temporal  variations  of  tropospheric 
winds,  a  42-hour  transhorizon  propagation  experiment  was  conducted  over 
a  164- km  asymmetrical  path  at  3.2  GHz.  The  propagation  path  extended 
from  Jackson  Butte,  California  (transmitter)  to  the  foothills  behind 
Stanford  University  (receiver).  A  rapid  data-gathering  system  (Cox  and 
Waterman,  1969)  used  in  previous  transhorizon  experiments  (Cox,  1967; 

Cox  and  Waterman,  1971;  Cianos,  1971;  Cianos  and  Waterman,  1973;  Water¬ 
man,  1973;  Cianos,  1978)  was  modified  to  remotely  probe  the  wind-field 
structure  in  the  common  volume  region  above  the  Livermore  Valley  about 
40  km  from  the  receiver. 

In  this  experiment  the  phase-coherent  system  employs  a  horizontal 
12-element  linear  receiving  antenna  array.  The  amplitudes  and  relative 
phases  of  individual  elements  are  sequentially  sampled  and  recorded  at 
a  40-Hz  rate.  An  effective  scan  period  of  5  ms  is  used  to  acquire  the 
data  from  all  twelve  elements.  Off-line  signal  processing  is  performed 
on  the  amplitude  and  phase  data  samples  to  "scan"  in  azimuth  a  0.3°  wide 
beam  through  a  3.88°  sector  and  to  obtain  corresponding  angle-of-arrival 
Doppler  spectra.  The  average  Doppler  shift  (spectrum  centroid)  is 
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computed  as  a  function  of  the  array  pointing  angle  and  then  interpreted 
in  terms  of  the  average  transverse  wind  speed  and  direction.  Results 
are  compared  with  in-situ  wind  observation  data  obtained  by  optically 
tracking  pilot  balloons  released  into  the  common  volume. 

The  experimental  data  presented  in  this  dissertation  were  obtained 
from  tests  conducted  on  April  21-22,  1978.  Ten  minutes  of  transhorizon 
radio  data  were  recorded  each  hour.  A  pilot  balloon  was  launched  hour¬ 
ly  between  0800  and  1800  hours  in  the  middle  of  every  10-minute  data 
recording  period. 


1 .4  Summary  of  Results 

Some  of  the  significant  results  of  this  research  are: 

(1)  Wind-Doppler  measurement  results  compare  most  favorably 
with  in-situ  PIBAL  data  in  lower  strata  of  the  common 
volume. 

(2)  Spectral  broadening  of  wind-produced  spectra  are  ob¬ 
served  during  afternoon  hours  when  the  atmosphere  is 
well  mixed. 

(3)  At  syirmetrical  off-axis  angles,  measured  upwind  and 
downwind  Doppler  spectra  are  observed  to  be  nearly 
antisymmetric  about  the  zero-Doppler  frequency.  This 
suggests  some  degree  of  spatial  uniformity  in  the 
transverse  winds. 

(4)  On-axis  the  Doppler  spectra  are  generally  symmetrical 
about  zero-Doppler  suggesting  the  absence  of  vertical 
wind  components  and/or  spatial  non-uniformities  in  the 
transverse  wind  velocity. 

(5)  Aircraft  trailing  vorticies  produce  spectral  broaden¬ 
ing  and  spatial  asymmetries  in  the  ambient  wind  pro¬ 
duced  spectra. 

(6)  Spectral  and  transit-time  analyses  of  aircraft-induced 
Doppler  permit  a  determination  of  the  spatial  extent, 
range,  and  velocity  of  trailing  vorticies. 
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(7)  Analytical  results  show  the  Doppler  shift  for  a 
scatterer  at  a  fixed  cross-path  position  to  be 
greater  on  a  longitudinally  asymmetrical  path  vis 
a  vis  a  symmetrical  path  of  equal  path  length. 
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Chapter  2 


THEORETICAL  DISCUSSION 

Signals  received  via  transhorizon  tropospheric  scattering  proces¬ 
ses  are  frequently  subject  to  rapid  fading  and  spectral  instabilities. 
This  is  a  consequence  of  the  interference  nature  of  the  Doppler-shifted 
components  associated  with  moving  off-path  scatterers.  Interference 
patterns  are  set  up  at  the  receiver  by  signal  components  scattered  off 
refractive  eddy  structures  transported  by  the  wind  through  regions  of 
the  atmosphere  mutually  illuminated  by  the  transmitting  and  receiving 
antenna  beams.  As  suggested  by  Atlas  et  al .  (1969),  atmospheric  struc¬ 
tures  on  the  upwind  side  of  the  beam  produce  positive  Doppler  frequency 
shifts  that  beat  with  negative-shifted  Doppler  components  associated 
with  eddies  on  the  downwind  side  of  the  beam.  These  Doppler  shifts 
progressively  increase  with  the  off-path  displacement  of  the  scattering 
region  and  have  been  shown  (Birkemeier  et  al.,  1968;  Atlas  et  al.,  1969) 
to  be  directly  proportional  to  the  cross-path  wind  speed  in  the  lower 
segments  of  the  common  volume.  A  single-scatterer  model  (Birkemeier 
et  al.,  1968)  can  be  used  to  illustrate  the  physical  concepts  involved 
and  to  show  the  influence  of  the  path  geometry  and  source  position  on 
the  basic  Doppler  relations. 

2.1  Single  Scatterer  Model 

In  the  analysis  presented  in  this  section  a  single  scatterer  is 
assumed  to  be  responsible  for  the  signal  observed  at  the  receiver.  As 
shown  in  Figure  2.1,  the  point  source  scatterer  is  located  at  the 


SCATTERING  SOURCE  TRANSMITTER 


Fig.  2.1.  SINGLE  SCATTERER  PATH  GEOMETRY. 


Intersection  of  two  perfectly  narrow  antenna  beams.  The  Cartesian  co¬ 
ordinate  system  used  in  this  model  has  its  origin  centrally  located 
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along  the  chordal  axis  between  the  two  radio  terminals.  The  x,  y,  and 
z  coordinates  define  respectively  the  longitudinal,  transverse,  and 
vertical  spatial  positions  of  the  scatterer.  The  total  ray-path  length, 

L  ,  from  the  transmitter  to  the  receiver  is  the  sum  of  the  respective 
path  lengths,  Lt  and  Lr  ,  from  the  transmitter  to  the  scatterer  and 
from  the  receiver  to  the  scatterer.  Denoting  the  chordal  distance  be¬ 
tween  the  terminals  as  2D  ,  the  total  ray-path  length  can  be  expres¬ 
sed  in  terms  of  the  position  coordinates  of  the  scatterer  as 

L  *  [(D+x)2  +  y2  +  z2]*5  +  [(D-x)2  +  y2  +  z2]*5  (2.1) 


If  we  assume  a  4/3  effective  earth  radius  (Bean  and  Dutton,  1966)  to  ac¬ 
count  for  the  refractive  bending  of  radio  rays  in  the  lower  atmosphere, 
we  may  use  straight  lines  for  the  ray  paths  in  applying  the  above  expres¬ 
sion. 

It  is  well  known  that  a  moving  source  produces  a  Doppler  shift, 
fj  ,  at  the  receiver  that  is  proportional  to  the  time  rate  of  change  of 
the  ray-path  length.  This  Doppler  principle  can  be  expressed  as 


1  dL 
X  dt 


(2.2) 


where  X  is  the  wavelength  of  the  transmitted  signal.  It  is  seen  that 
if  the  effective  path  length  decreases  (increases)  with  time,  a  positive 
(negative)  Doppler  shift  will  be  observed  at  the  receiver.  The  Doppler 
shift  can  also  be  interpreted  as  the  rate  at  which  the  scatterer  moves 
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across  surfaces  of  constant  path  length.  Contours  of  constant  path 
length  are  ellipsoidal  surfaces  (Birkemeier  et  al . ,  1968)  with  the  radio 
terminals  as  foci.  If  the  scatterer  moves  with  a  velocity,  7s  ,  an 
equivalent  Doppler  expression  can  be  written  as 

fd  ■  -K  • 7L  (2-3» 


Upon  substituting  the  expression  given  in  (2.1)  for  the  total  ray-path 
length,  the  Doppler  shift  from  a  single  moving  scatterer  becomes  expres 
sed  as 


[ux  +  vy  +  wz]L  +  uD[(y2-2xD)is  -  (y2+2xD)'1] 
(y2+2xD)^  (Y2-2xD)Js 


(2.4) 


where  In  the  above  y2  =  [x2  +  y2  +  z2  +  D2]  and  u,  v,  and  w  are  re¬ 
spectively  the  longitudinal,  transverse,  and  vertical  velocity  components 
of  the  moving  scatterer.  From  the  general  result  in  (2.4)  specific  Dop¬ 
pler  expressions  for  symmetrical  and  asymmetrical  paths  will  be  derived 
to  illustrate  the  influence  of  the  scatterer' s  longitudinal  position. 


2.2  Symmetrical  Path  Case 

On  longitudinally  symmetrical  paths,  the  scattering  regions  pri¬ 
marily  responsible  for  the  received  signal  are  located  in  close  proxi¬ 
mity  of  the  mid-path  plane.  As  a  result,  a  moving  point-source  scat¬ 
terer  has  relatively  small  longitudinal  displacements,  i.e.  x  «  D  and 
is  generally  confined  to  scattering  regions  in  the  common  volume  where 
D2  »  [x2  +  y2  +  z2]  .  Applying  a  first-order  approximation  to  the  terms 
(y2+2xD)J5  and  (y2-2xD)?s  in  (2.4),  the  Doppler  shift  produced  by  a 
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moving  scatterer  on  a  symmetrical  path  becomes  expressible  as 

(2.5) 

This  result  arrived  at  by  Birkemeier  et  al .  (1968)  has  been  applicable 
to  previously  reported  transhorizon  wind-sensing  experiments  since  such 
experiments  have  heretofore  been  conducted  on  near-symmetrical  paths. 

Equation  (2.5)  can  be  further  simplified  by  considering  the  geome¬ 
try  of  the  symmetrical  path  profile  in  Figure  2.2  and  the  relative  Lon- 
tributions  of  the  wind  velocity  components.  It  can  be  shown  that  the 

term  (1-4D2/L2)  is  small  and  equivalent  to  4z2/L2  =  4h2/D2  where  z 

mm  m 

and  h^  are  respectively  the  mid-path  heights  of  the  scatterer  above  the 
chordal  axis  and  the  earth's  surface.  Since  the  longitudinal  displace¬ 
ment,  x  ,  of  the  scatterer  is  also  small,  the  Doppler  contribution  from 
the  longitudinal  wind  can  be  neglected,  the  longitudinal  motion  produ¬ 
ces  negligible  Doppler  effects  because  of  the  elongated  shapes  of  the 
ellipsoidal  surfaces  of  constant  path  length.  A  scatterer  moving  in 
the  longitudinal  direction  cuts  across  the  surfaces  at  a  much  slower 
rate  than  one  moving  in  the  transverse  or  vertical  directions.  The  ef¬ 
fect  of  the  vertical  wind  can  in  most  instances  be  neglected  since  its 
magnitude  is  normally  much  less  than  that  of  the  horizontal  wind.  It 
then  becomes  evident  that  the  Doppler  shift  on  a  symmetrical  path  is 
linearly  proportional  to  the  cross-path  displacement,  y  ,  of  the  mov¬ 
ing  scatterer  and  may  be  explicitly  given  as 


_ 4  ux(l-4D2/L2)  +  vy  +  wz 

XL  1  -  16x2D2/LH 
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The  Doppler  shift  may  also  be  expressed  in  terms  of  the  azimuthal  angu¬ 
lar  position,  a  ,  of  the  scatterer.  The  Doppler  relationship  is  then 
written  as 

fd  ’  ■  rsin(a)  (2<7) 


which,  for  small  angles,  can  be  approximated  with  negligible  error  as 


2yg 

X 


(2.8) 


2.3  Asymmetrical  Path  Case 

Consider  now  the  path  profile  shown  in  Figure  2.3  where  the  prin¬ 
cipal  scattering  region  in  the  common  volume  is  skewed  significantly  to 
one  side  of  the  mid-path  plane.  We  begin  this  treatise  by  again  looking 
at  the  terms  (y2+2xD)^  and  (y2-2xD)^  in  the  general  Doppler  shift  ex¬ 
pression  given  in  (2.4)  but  view  them  from  the  standpoint  of  an  asym¬ 
metrical  path.  For  this  situation,  it  is  assumed  that  the  longitudinal 
displacement,  x  ,  is  such  that  (D-x)2  »  y2  +  z2  for  all  values  of 
x  ,  y  ,  and  z  within  the  common  volume.  With  this  constraint,  a  first- 
order  approximation  to  the  terms,  (y2+2xD)^  and  (y2-2xD)^  ,  results  in 
a  Doppler  shift  on  an  asymmetrical  path  that  can  be  expressed  as 


L 

XD2 


2ux(l -2D/L)  +  vy  +  wz 


1  -  x2/D2 


(2.9) 


This  result  reduces  to  that  for  the  symmetrical  path  case  when  x  «  D  . 
Hence,  the  Doppler  shift  expression  given  in  (2.9)  is  a  more  general  re¬ 
sult  that  is  applicable  to  both  symmetric  and  asymmetric  path  geometries. 
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A  simplification  of  (2.9)  can  be  made  by  introducing  a  path  asym¬ 
metry  factor  resulting  from  the  terrain  profile  along  the  path  and  the 
relative  heights  of  the  radio  terminals.  The  asymmetry  factor,  s  ,  is 
defined  as  the  ratio  of  the  longitudinal  displacement  of  the  scattering 
source  to  the  chordal  distance  from  the  mid-path  plane  to  either  of  the 
radio  terminals,  i.e.  s  =  x/D  .  It  Is  zero  for  a  symmetrical  path. 

When  the  scattering  region  is  skewed  to  one  side  of  the  mid-path  plane, 
however,  the  magnitude  of  the  asymmetry  factor  becomes  greater  than  zero 
but  always  remains  less  than  unity.  Using  this  factor,  the  Doppler  shift 
becomes 

2us(l-2D/L)  +  vy/D  +  wz/D  ^ 

1  -  s2 


The  term  (1-2D/L)  is  approximately  equivalent  to  (y2+z2)/2(D2-x2)  . 
Thus,  using  the  geometry  in  Figure  2.1,  the  Doppler  shift  can  then  be 
expressed  in  terms  of  the  azimuthal,  a  ,  and  vertical,  3  ,  angular 
positions  of  the  scatterer  as 


f  =  - 

"d  XD 


us(arat  +  Br8t)  var  +  wBr 


1  -  s' 


1  +  s 


(2.11) 


where  small  angle  approximations  with  respect  to  the  receiver  and  trans- 
mitter  ere  used,  i.e.  ^  =  gfl^y  ,  ot  «  gfl^y  ,  6r  -  , 

and  -  ~p^-+sy  .  For  equal  transverse  and  longitudinal  wind  veloci¬ 
ties,  the  longitudinal  Doppler  contribution  again  remains  relatively 
small  and  can  be  neglected.  If  the  vertical  wind  contribution  is  also 
neglected  and  the  first-order  approximation,  L  =  2D  ,  is  substituted 
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for  the  ray-path  length,  the  Doppler  shift  expression  (2.10)  simplifies 
to 


>  „  .  2vy 
d  XDO-s2} 


(2.12) 


which  can  be  rewritten  in  terms  of  the  azimuthal  angular  offset,  ar  , 
from  the  receiver  as 

2va 

fd  ■  -  xrr^y  <2-131 

Comparing  the  above  results  with  those  for  the  symmetrical  path  case, 
equations  (2.6)  and  (2.8),  it  is  seen  that  the  Doppler  shift  for  a  scat- 
terer  at  a  given  cross-path  displacement  is  greater  on  an  asymmetrical 
path  than  on  a  symmetrical  path.  This  is  consistent  with  the  interpreta¬ 
tion  of  the  Doppler  shift  as  the  rate  at  which  a  scatterer  cuts  across 
surfaces  of  constant  path  length.  The  separation  between  these  ellipsoi¬ 
dal  surfaces  is  greatest  in  the  mid-path  plane  and  becomes  smaller  to¬ 
wards  either  of  the  two  terminals.  Consequently,  a  scatterer  moving 
through  the  common  volume  region  on  an  asymmetrical  path  cuts  across  the 
equiphase  surfaces  at  a  more  rapid  rate  than  on  a  symmetrical  path. 

On  the  asymmetrical  Jackson-Stanford  path  used  in  the  experiment 
discussed  in  this  dissertation,  the  principal  scattering  region  is  in  the 
vicinity  of  a  transverse  plane  about  40  kilometers  from  the  receiver. 

The  average  longitudinal  displacement,  x  ,  of  a  scatterer  on  this  path 
is  such  that  s  -  %  .  The  analytical  expression  for  the  Doppler  shift 
applicable  to  this  path  is 
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(2.14) 


4yar 

3X 


where  a  Is  measured  in  radians.  If  the  angle  measurement  is  conver- 
r 

ted  to  degrees  and  the  numerical  value  (9.4  x  10"2m)  is  substituted  for 
X  ,  a  more  simplified  expression  for  the  Doppler  shift  on  the  Jackson- 
Stanford  path  can  be  written  as 

va 


fd  =  ‘ 


(2.15) 


where  in  the  above  expression  v  is  units  of  m/s. 

2.4  Multiple  Scattering  Effects 

In  the  previous  sections,  a  single  scattering  source  was  assumed  and 
as  a  result,  Doppler  effects  were  limited  to  single  frequency  shifts.  A 
multiplicity  of  scatterers  randomly  distributed  throughout  an  illuminated 
region  in  the  atmosphere  produces  a  spectrum  of  Doppler  shifts.  The 
spectral  composition  of  the  Doppler  shifts  is  dependent  upon  the  spatial 
distribution  and  the  relative  velocities  of  the  individual  scatterers. 
Each  scatterer  gives  rise  to  a  discrete  Doppler  shift  in  accordance  with 
the  analysis  developed  earlier.  A  number  of  the  scatterers,  however,  may 
contribute  to  the  same  spectral  component.  The  extent  of  the  contribu¬ 
tions  from  the  many  off-path  scatterers  is  reflected  in  the  shape  of  the 
Doppler  spectrum,  which  is  dependent  upon  the  relative  strength  of  the 
received  signal  as  a  function  of  the  scatterers'  cross-path  positions. 

The  signal  power  received  from  a  particular  azimuthal  direction  is  di¬ 
rectly  proportional  to  the  product  of  the  net  antenna  gain  and  the  re¬ 
flectivity  or  scatter  cross-section  of  the  medium.  Thus,  the  Doppler 
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spectrum  arising  from  a  distribution  of  scatterers  is  dependent  upon  the 
antenna  beam  configuration  and  the  basic  propagation  mechanism. 

An  average  Doppler  shift,  f^  ,  can  be  attributed  to  a  mean  flow 
or  drift  of  the  atmospheric  scatterers.  Birkemeier  et  al .  (1968)  have 
specified  analytically  the  average  Doppler  shift  for  antennas  synchro¬ 
nously  aligned  in  azimuth  at  an  angle,  a  ,  as 

d 

fj  *  -  Y~J'aG[a-aa)U[a)dia  / J' I G(a-aa )W(a)da  (2.16) 

where  G(a-a  )  is  the  combined  antenna  gain  function  and  W(a)  is  the 
scattered  power  per  unit  azimuthal  angle,  this  result,  however,  applies 
to  a  symmetrical  path.  A  more  general  result  incorporates  the  path  asym¬ 
metry  factor,  s  ,  and  can  be  given  in  terms  of  the  angular  offset, 

a  ,  measured  from  the  receiver  location,  as 
r 

fd  =  x(l+s]  /*(a-ar)W(a)da  /J G(a-ar )W(a)da  (2.17) 

In  these  expressions  for  the  average  Doppler  shift,  various  atmospheric 
models  can  be  applied  to  account  for  the  angular-dependence  scattering 
function,  W(a)  .  The  classical  isotropic  turbulence  model  (Tatarski, 
1961)  predicts  a  minus  11/3  power-law  variation  in  the  scattering  cross- 
section  as  a  function  of  scattering  angle,  9  .  The  scattering  angle  is 

defined  as  the  angle  between  the  transmitted  and  received  ray-path  di¬ 
rections.  For  the  isotropic  turbulence  model,  the  scattering  cross- 

_  li 

section  is  proportional  to  sin  3  (0/2)  .  Birkemeier  et  al .  (1977) 

have  shown  that  for  small  scattering  angles  near  the  mid- path  plane 
9  =  2>lo}  +  02  .  This  expression  for  the  scattering  angle  is  valid  for 
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near-symmetric  propagation  paths.  However,  a  more  general  result  that 
Is  applicable  to  asymmetrical  paths  as  well  shows  the  scattering  angle 
to  be  given  approximately  as 

6  -  +  3^  +  Vat  +  &t  (2.18) 

Using  the  small  angle  approximations  developed  earlier,  the  above  expres¬ 
sion  can  be  written  as  a  function  of  the  path  asymmetry  factor,  s  ,  in 
terms  of  the  angular  displacements,  and  8^  ,  from  the  receiver  as 

8  *  t4t  V“!r  ♦  (2-'9) 

An  equivalent  expression  for  the  Jackson-Stanford  path,  which  has  a  mini¬ 
mum  scattering  angle  of  0.84®  occurring  at  an  elevation  angle  of  0.25°, 
is 

e  =  0.5  +  J  yjai*  +  ej  (2.20) 


where  8  >  0.25®  . 

r  — 

In  addition  to  the  azimuthal  variations,  the  above  relationships 
show  the  extent  of  the  influence  of  the  vertical  positions  of  the  scat¬ 
tered  on  the  scattering  function.  The  average  Doppler  shift  from  a 
given  azimuthal  direction  is  a  spatially  averaged  result  of  the  trans¬ 
verse  motions  of  scatterers  at  various  altitudes  within  the  common 
volume.  The  wind-produced  drift  of  the  scatterers  generally  increases 
with  altitude.  However,  the  turbulence  model  assumed  produces  a  rapid 
decrease  in  the  scattered  power  with  elevation  angle.  This  suggests 
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the  average  Doppler  shift  is  primarily  associated  with  the  motion  of 
scatterers  In  the  lower  segments  of  the  illuminated  scattering  region. 
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Chapter  3 


DESCRIPTION  OF  EXPERIMENT 

3.1  Propagation  Path 

The  experiment  described  in  this  report  was  conducted  over  the  pro¬ 
pagation  path  shown  in  Figure  3.1.  The  164-km  (1C?  miles)  path  is  align¬ 
ed  along  a  SW-NE  (230°  -  50°)  azimuthal  bearing  between  Jackson  Butte 
(transmitter  site)  and  Stanford,  California  (receiver  site).  The  trans¬ 
mitter  was  situated  on  a  725-meter  (2375  ft)  peak  just  east  of  Jackson, 
California  and  the  receiver  was  located  110  meters  (360  ft)  above  sea  le¬ 
vel  in  the  foothills  behind  Stanford  University.  A  topographical  profile 


Fig.  3.1.  PROPAGATION  PATH  AND  SURROUNDING  AREA. 
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Fig.  3.2.  TOPOGRAPHICAL  PROFILE  OF  PROPAGATION  PATH. 


(Figure  3.2)  of  the  path  shows  the  existence  of  three  mountain  ridges 
that  are  nearly  perpendicular  to  the  path  between  the  two  terminals. 
These  ridges  determine  the  local  radio  horizons  of  the  transmitting  and 
receiving  antenna  beams  which  are  aligned  so  as  to  mutually  illuminate  a 
region  from  about  425  meters  (1400  ft)  to  approximately  2  km  (6600  ft) 
above  sea  level.  The  path  geometry  introduces  a  minimum  scattering  an¬ 
gle  of  0.84° 

3.2  Instrumentation 

3.2.1  Data-Gathering  Array 

Angle-of-arrival  Doppler  measurements  are  derived  from  the  indivi¬ 
dual  amplitude  and  relative  phase  sampled  outputs  of  a  12-element  data 
gathering  antenna  array.  The  array,  which  operates  at  3.2  GHz,  is  shown 
in  Figure  3.3  and,  except  for  the  horizontal  orientation  and  polariza¬ 
tion,  is  identical  to  the  vertical  array  used  in  previous  experiments  at 
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HEIGHT  (km) 


Stanford  (Cox,  1967;  Cox  and  Waterman,  1969  and  1971;  Cianos,  1971; 

Cianos  and  Waterman,  1973;  Waterman,  1973;  and  Cianos,  1978).  The  data- 
gathering  array  is  aligned  nearly  broadside  to  the  propagation  path  and 
is  structurally  supported  beneath  the  2nd,  6th,  and  10th  elements.  Each 
element  is  horizontally  polarized  and  has  a  5°  half-power  beamwidth  which 
determines  the  vertical  extent  of  the  array  fan  beam.  When  uniformly  il¬ 
luminated,  the  array  produces  in  azimuth  a  0.3®  wide  beam  which  can  be 
scanned  through  a  3.88°  azimuthal  sector.  A  summary  of  the  physical  and 
operational  characteristics  of  the  array  is  outlined  in  Table  3.1. 


Table  3.1 

SUMMARY  CHARACTERISTICS  OF  DATA-GATHERING  ARRAY 


Physical 

Characteristics 

Number  of  elements 

Element  spacing 

Element  diameter 

Effective  aperture  length 

Twel ve 

1.4  m  (14.8  wavelengths) 

1.2  m  (13.0  wavelengths) 

15.2  m  (162.5  wavelengths) 

Operating  frequency 

3.2  GHz 

Operational 

Polarization 

Horizontal 

Characteristics 

Half-power  beamwidth* 

0.3°  AZ/5®  EL 

Grating  lobe  separation* 

3.88° 

Sidelobe  levels* 

<  -  13  dB 

*  for  a  uniformly  illuminated  conventional  additive  array 
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3.2.2  Array-Receiver  Channels 


Twelve  separate  receiving  channels  connect  to  the  twelve  separate 
elements  in  the  array.  A  typical  receiving  channel  is  illustrated  in 
Figure  3.4.  Each  channel  receives  simultaneously  the  transhorizon  sig¬ 
nal  at  3.2001  GHz  and  a  stable  1 ine-of-sight  reference  signal  at  3.2000 
GHz.  A  3.2300-GHz  local  oscillator  signal  (common  to  all  12  channels) 
is  mixed  with  both  the  transhorizon  and  reference  signals  in  mixer- 
preamps  mounted  directly  behind  each  array  element.  The  down-converted 
transhorizon  (29.9-MHz)  and  reference  (30.0-MHz)  signals  are  separated 
by  2-kHz  wide  crystal  filters.  The  30-MHz  IF  reference  is  amplified 
and  then  used  as  a  second  local  oscillator  which  mixes  with  the  29.9-MHz 
transhorizon  IF  signal.  The  resultant  100-kHz  second  IF  signal  is  nar¬ 
row-band  filtered  with  a  23-Hz  wide  crystal  filter  which  effectively 
sets  the  noise  bandwidth  of  the  receiving  channel.  After  further  ampli¬ 
fication  the  narrow-band  100-kHz  signal  is  passed  through  a  variable 
phase  shifter  and  a  variable  gain  amplifier  which  are  adjusted  and  set 
during  calibration.  This  signal  is  then  continuously  available  from 
each  array-element  receiving  channel. 

Sampling  of  the  twelve  100-kHz  signals  is  performed  by  a  computer- 
controlled  electronic  commutator.  A  block  diagram  of  the  amplitude  and 
phase  sampling  equipment  is  shown  in  Figure  3.5.  Every  25  ms  the  com¬ 
mutator  is  sequentially  stepped  through  all  12-element  positions.  The 
commutator  takes  five  milliseconds  to  completely  scan  all  twelve  ele¬ 
ments.  The  commutator  output  is  bandpass  filtered  to  reduce  switching 
transients.  The  sample  100-kHz  output  is  phase  compared  with  a  stable 
100-kHz  reference  standard.  Additionally,  the  signal  is  applied  to  two 
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Fig.  3.4.  BLOCK  DIAGRAM  OF  AN  ARRAY  ELEMENT  RECEIVING  CHANNEL. 

Parameters 


1st  Mixers  (part  of  LEL  type 
SBC-6 

Noise  figure:  9  to  10  dB 
(depends  on  unit) 

30-MHz  IF  Preamplifiers  (part 
of  LEL  type  SBC-6) 

Gain:  25  to  30  dB  RF  to  IF 
(depends  on  unit) 

30-MHz  IF  Amplifiers  (LEL  type 

iron 

Adjustable  gain:  20  to  70  dB 
Bandwidth:  3  MHz 

29.9-MHz  and  30,0-MHz  Filters 
(Hughes  special  design) 

Bandwidth:  2  kHz 
Quartz  crystals 

2nd  Mixers 

Balanced  diodes  for  LO  (refer¬ 
ence  signal) 

Noise  Suppression 


100-kHz  Filters  (Blackhawk  Networks 
Corporation) 

Bandwidth:  23  Hz 
Quartz  Crystals 

Phase  Shifters  (Nil sen  V42-01  Vari- 
ogon) 

Adjustable  0°  to  360°  continuously 
Loss  independent  of  phase  setting 

100-kHz  Output  Amplifiers 


Low  output  impedance:  5  to  10  ft 
Phase  shift  independent  of  gain 
setting 

Gain  adjustable  0  to  20  dB 


All  active  circuit  elements  are 
solid-state  devices 
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Fig.  3.5.  BLOCK  DIAGRAM  OF  AMPLITUDE  AND  PHASE-SAMPLING  EQUIPMENT. 


Parameters 


Electronic  Commutator 

2-pole  12-position  switch 
Scans  array  at  100  scans  (over 
12  elements)/sec 
Settling  time:  «  0.3  msec 

100-kHz  Bandpass  Filters 
Bandwidth:  30  kHz 

Phase-Angle  Computer  (Wiltron 

Mode  13  51) 

Measures  phase  to  1°  or  better 
over  60  dB  signal  variation 
in  either  or  both  channels 

Lowpass  Filters 

Active  2-pole  filter  using  RC 
feedback  around  an  opera¬ 
tional  amplifier 
RC  feedback  set  to  optimize 
rise  time  at  switch  rate 
used 

dc  offset:  <  1  mV 


Linear  Amplitude  Detector 

Linear  to  1%  of  full  scale  over 
dynamic  range  >  40  dB 


dc  Amplifiers 

dc  offset:  <  1  mV 
Feedback  around  operational 
amplifier  for  good  linearity 


Frequency  Standard  (James  Kniqht 
FS-llOOT 


100  kHz  output 

Aging  rate  of  1  to  2  parts  in 
1010/day  (1  to  2-1 0“5  Hz/da y 
at  100  kHz) 
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linear  amplitude  detectors  one  of  which  is  used  exclusively  for  ampli¬ 
tude  calibration.  The  detected  amplitude  and  relative  phase  data  samples 
are  low-pass  filtered  and  made  available  for  continuous  monitoring  on  an 
oscilloscope  before  being  processed  for  data  recording. 

A  block  diagram  of  the  data  collection  equipment  is  shown  in  Figure 
3.6.  The  detected  amplitude  and  phase  data  samples  are  multiplexed  and 
converted  to  12-bit  digital  words  under  the  control  of  a  general  purpose 
PDP-8  computer.  The  12-bit  digital  format  allows  for  a  maximum  of  4096 
quantization  levels  in  the  analog-to-digital  (A/D)  conversion  process. 


advance  pulse 


Fig.  3.6.  BLOCK  DIAGRAM  OF  DATA-COLLECTION  EQUIPMENT. 


Parameters 


A/D  Converter  and  Multiplexer 
(DEC  AFOT) 

Multiplexer  switching  speed: 
<  2  ysec 

A-D  word  length:  12  bits 
(3096  levels) 

35  psec  conversion  time 
+0.25%  max  switching  point 
error  for  +1/2  least  sig¬ 
nificant  bit  of  quantiza¬ 
tion  error 


Computer  (DEC  PDP-8)  with  Hardware 
Arithmetic  Unit 

1.5  psec  memory  cycle  time 
4096  word  fast  memory 
12-bit  word  length 

Magnetic  Tape  Units 

DEC  TU-20  (IBM  Compatible) 

DEC  TU-55  (Program  DECtape) 


DUMMY 

LOAD 


Fig.  3.7.  BLOCK  DIAGRAM  OF  TRANSMITTER 
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The  buffer  amplifiers  preceding  the  A/D  converter  are  adjusted  so  that 
the  full  dynamic  range  of  the  converter  is  utilized.  In  converting  the 
phase  measurements,  a  range  of  -  180®  to  +  180®  corresponds  respectively 
to  converter  readouts  of  0  to  12  bits  (0  to  4095  levels);  in  converting 
the  amplitude  measurements,  the  gains  are  set  so  that  no  input  signal 
corresponds  to  0  bits  (0  level)  and  a  fixed  calibration  signal  corre¬ 
sponds  to  11  bits  (2048  level).  After  analog-to-digital  conversion,  the 
amplitude  and  phase  data  samples  are  written  and  stored  on  magnetic  tape. 

3.2.3  Transmitter  Assembly 

A  block  diagram  of  the  transmitting  system  is  shown  in  Figure  3.7. 
The  transmitter  is  essentially  the  same  as  that  used  previously  by  Cox 
(1967),  Cox  and  Waterman  (1971),  Cianos  (1971  and  1978),  Cianos  and 
Waterman  (1973),  and  Waterman  (1973).  The  transportable  transmitting 
antenna  shown  in  Figure  3.8  uses  the  same  2.4-m  (8-ft)  parabolic  reflec¬ 
tor  as  before  but  employs  a  horizontally  polarized  feed.  The  transmit¬ 
ter  has  available  a  maximum  output  power  of  20  watts;  however,  during 
the  actual  tests,  the  nominal  operation  for  the  transmitter  was  between 
10  and  12  watts.  Some  modifications  were  made  in  the  output  power¬ 
monitoring  assembly  so  that  power  level  adjustments  could  be  performed 
while  the  transmitter  was  in  a  non-radiating  mode.  These  modifications 
included  the  use  of  a  dummy  load  and  additional  coaxial  switches  to 
permit  the  transmitter  signal  to  be  monitored  under  loaded  non-radi¬ 
ating  conditions  during  calibration  and  alignment  periods. 


Fig.  3.8.  TRANSMITTING  ANTENNA  AND  EQUIPMENT  VAN. 


3.3  Experimental  Procedures 


3.3.1  Data  Collection 

Data  was  collected  and  recorded  during  a  continuous  42-hour  experi¬ 
ment  conducted  between  0000  hours  on  21  April  1978  and  1830  hours  on  22 
April  1978.  Ten  minutes  of  transhorizon  data  were  recorded  every  hour 
beginning  approximately  on  the  hour.  The  sampled  amplitude  and  phase 
data  from  each  array  element  were  recorded  in  12-bit  digital  format  on 
magnetic  tape  at  a  40-Hz  rate.  A  five  millisecond  scan  period  was  used 
to  sample  and  record  the  amplitude  and  phase  data  from  all  twelve  ele¬ 
ments.  The  data  samples  were  recorded  In  the  sequence  a  ,4>  ;  a  ,<J>  ;  .  . 

11  2  2 

.  .  a  *4*  • 

12  12 

A  calibration  signal  transmitted  from  a  1 ine-of-sight  source  about 
a  mile  from  the  receiver  along  the  path  was  recorded  for  1.6  seconds  (64 
scans)  before  every  ten-minute  data  recording  period.  Immediately  fol¬ 
lowing  the  recording  of  transhorizon  data,  the  transmitter  was  turned  off 
and  the  calibration  signal  again  turned  on.  The  amplitude  and  phase  cal¬ 
ibration  levels  of  each  array  element  were  then  monitored  with  the  digi¬ 
tal  voltmeter  (Figure  3.5)  and  read  into  the  log  while  the  commutator  was 
manually  stepped  through  all  twelve  element  positions. 

During  actual  data  recording,  the  sampled  amplitudes  and  phases  of 
all  twelve  elements  were  continuously  monitored  on  an  oscilloscope  (Fig¬ 
ure  3.5).  Any  signal  anomalies  that  could  be  attributed  to  aircraft. 
Interference,  etc.  were  noted  in  the  log.  Visual  sightings  of  aircraft 
In  the  vicinity  of  the  common  volume  were  logged  at  the  transmitter, 
receiver,  and  PIBAL  launch  sites. 
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3.3.2  System  Calibration 

Prior  to  each  transhorizon  data  recording  period,  the  transmitter 
was  completely  checked  and  aligned;  and  all  array  element  receiving  chan¬ 
nels  were  calibrated.  The  transmitter  frequency  standard  (Figure  3.5) 
was  adjusted  each  day  during  the  quiet  early  morning  hours  and  set  to 
within  2  parts  in  1011  (0.064  Hz  at  3.2  GHz)  of  the  frequency  standard 
used  to  generate  the  reference  and  calibration  signals.  In  the  cali¬ 
bration  and  alignment  of  the  individual  receiving  channels,  a  number  of 
critical  checks  were  made  to  ensure  proper  operation  during  the  data- 
gathering  portions  of  the  experiment.  The  mixer  currents,  the  amplifier 
gains,  and  the  phase-lock  loop  status  of  the  receiving  system  were  rou¬ 
tinely  checked  and  if  necessary,  adjusted  before  each  data  recording 
period. 

During  calibration  periods,  the  transhorizon  signal  was  turned  off 
and  a  3.2001 -GHz  calibration  signal  was  transmitted  from  a  stable  source 
about  1.6  kilometers  from  the  receiver.  This  signal  was  remotely  con¬ 
trolled  from  the  receiving  site  and  was  used  to  properly  set  the  vari¬ 
able  100-kHz  amplifiers  and  phase-shifters  in  each  receiving  channel. 

With  the  calibration  signal  on,  the  commutator  was  manually  advanced 
through  all  twelve  element  positions  and  the  detected  amplitude  and 
adjacent  element  phase-difference  signals  were  monitored  with  the  digi¬ 
tal  voltmeter.  The  variable  phase-shifters  and  amplifiers  were  then 
adjusted  on  an  element  by  element  basis  for  predetermined  calibration 
readings  on  the  digital  voltmeter. 

The  predetermined  phase-difference  and  amplitude  calibration  set- 
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tings  were  obtained  from  the  results  of  a  series  of  1  l'ne-of-slght  (LOS) 
tests  In  which  the  transmitter  was  placed  at  selected  azimuthal  posi¬ 
tions  about  15  kilometers  from  the  receiving  array.  With  the  calibration 
signal  off  and  the  transmitter  on  (at  reduced  power),  the  individual  am¬ 
plitudes  were  set  uniformly  across  the  array  by  adjusting  the  gains  of 
the  100-kHz  amplifiers.  Adjacent  element  phase-differences  were  set  with 
the  variable  phase-shifters  to  values  calculated  from  the  signal's  known 
angle-of-arrival .  After  the  phase  and  amplitude  adjustments  were  made, 
the  transmitter  was  turned  off  and  the  LOS  calibration  source  turned  on. 
Phase  and  amplitude  readings  were  then  obtained  for  all  twelve  elements 
using  the  digital  voltmeter. 

The  above  procedures  were  repeated  for  three  different  1 ine-of-sight 
transmitter  positions.  Results  from  two  of  the  azimuthal  positions  were 
averaged  and  used  to  obtain  the  final  amplitude  and  phase-difference  cal¬ 
ibration  settings.  The  maximum  variation  in  the  adjacent  element  phase- 
differences  observed  between  the  two  locations  was  14°  and  the  rms  vari¬ 
ation  was  7°.  In  the  amplitude  readings,  the  maximum  variation  was  2  dB 
and  the  rms  variation  was  0.7  dB. 

3.3.3  PIBAL  Wind  Monitoring 

Pilot  balloons  released  into  the  common  volume  from  a  launch  site 
near  Livermore  were  optically  tracked  during  daylight  hours  with  a  single 
theodolite  system  to  obtain  vertical  profiles  of  the  horizontal  wind- 
field.  A  pilot  balloon  (PIBAL)  was  launched  every  hour  at  five  minutes 
past  the  hour  between  0805  and  1805  on  each  of  the  two  data  recording 
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days.  The  balloons  were  released  by  members  of  the  San  Jose  State  Uni- 
verstty  Meteorology  Department  from  a  position  beneath  the  common  volume 
about  0.9®  SE  of  the  great-circle  bearing  and  165  meters  (540  ft)  above 
mean  sea  level.  The  PIBAL  ascent  rate  was  approximately  200  meters  (650 
ft)  per  minute.  Azimuth  and  elevation  position  readings  were  taken  every 
30  seconds  and  each  PIBAL  was  tracked  to  an  altitude  of  approximately 
2.4  km  (7870  ft)  above  sea  level.  An  algorithm  used  by  the  meteorology 
group  reduced  the  PIBAL  azimuth  and  elevation  information  to  horizontal 
wind  velocity  profiles. 
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Chapter  4 


ANALYTICAL  METHODS 

Two  basic  signal  processing  steps  are  used  to  obtain  the  angle-of- 
arrival  Doppler  spectra  from  measured  amplitude  and  phase  data  samples. 
First,  a  digital  Fourier  transform  is  performed  in  the  time  domain  on  the 
sampled  output  signals  from  the  individual  array  elements.  The  ampli¬ 
tudes  and  phases  of  the  resultant  spectra  are  preserved  and  used  to 
characterize  the  array  element  output  signals  in  the  frequency  domain. 

A  spatial  transform  is  then  performed  in  the  array  plane  on  the  element 
spectra  to  obtain  an  angular  Doppler  spectrum  for  the  array.  In  this 
process,  the  individual  element  spectra  are  combined  using  conventional 
additive  array  techniques  to  maximize  the  output  spectra  for  a  particu¬ 
lar  angle-of-arrival .  An  average  Doppler  frequency  shift  (spectrum  cen¬ 
troid)  is  computed  as  a  function  of  the  angular  "pointing"  direction  of 
the  array  with  respect  to  the  great-circle  (zero-Doppler)  plane.  The 
apparent  source  direction  together  with  the  computed  average  Doppler 
shift  is  used  to  infer  an  average  transverse  wind  velocity. 

4.1  Single  Element  Spectra 

The  measured  element  amplitude,  an(k)  ,  and  relative  phase, 

4>n (k)  ,  data  samples  are  used  to  construct  the  single  element  Doppler 

spectra.  Each  spectral  component  is  evaluated  from  a  set  of  K  con¬ 
tiguous  data  samples.  The  in-phase  and  out-of-phase  quadrature  com¬ 
ponents  of  the  element  spectra  are  computed  in  accordance  with  the 
following 
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where :  n  =  1 ,  2 , 

....  12 

The  amplitudes,  A~(f.)  ,  and  phases,  Y  „(f.)  ,  of  the  spectral  com- 

nN  I  n  i\  I 

ponents  are  computed  from  the  quadrature  phasor  components 

AnK(fi )  =  V  sXK(n*fi )  +  SYK(n»fi }  (4*3) 

fnK(f1)  =  tan^rSy^n.f^/^^n.f.)]  (4.4) 

These  parameters  are  obtained  for  100  Doppler  frequencies  in  a  preselect¬ 
ed  spectral  observation  window.  The  spectral  window,  which  is  less  than 
the  receiver  bandwidth  (23  Hz),  is  centered  about  zero  Doppler  to  permit 
simultaneous  observation  of  the  upwind  (positive)  and  downwind  (negative) 
Doppler  shifts. 


4.2  Array  Angular  Spectra 

The  array  angular  spectra  are  derived  from  the  amplitudes,  A  Jf,. )  , 

n  i\  l 

and  phases,  ¥n K(f ^ )  .  of  the  individual  element  spectra.  A  progressive 

linear  phase-shift  across  the  array  is  applied  to  the  spectral  components 
of  the  individual  element  spectra.  Each  component  is  phase-shifted  in 
accordance  with  conventional  additive  array  techniques.  The  quadrature 
phasor  components  of  the  resultant  angular  spectra  are  then  computed  for 
selected  beam-pointing  angles,  ,  between  +1.94°  as  follows 


40 


SXK(Vf1> 

SYK*°r’V 


Wcos  Kk(V  *  ("-’>1^]  t4.5) 

ZJflnK(f1)s1"  [W  ’  <n-'  K  1  <«S> 

where:  y  =  phase  shift  constant 

=  92.71  phase  degrees/beam  angle  degree 


A  positive  (negative)  angular  direction  corresponds  to  a  SE  (NW)  angle- 
of-arrival  with  respect  to  the  great-circle  bearing.  The  Doppler  or  po¬ 
wer  spectra,  P(otr » f ^ )  »  are  computed  for  a  particular  angular  direction 
by  summing  the  squares  of  the  quadrature  phasor  components,  i.e. 

P(ar,f.)  =  SjjK(a,,f.)  +  S5K(ar,f.)  (4.7) 


These  angular  spectra  are  essentially  unweighted  since  each  spectral  com¬ 
ponent  is  evaluated  from  the  same  set  of  element  amplitudes  irrespective 
of  its  source  direction.  The  Doppler  spectra,  however,  can  be  weighted 
to  account  for  the  angl e-of-arrival  variations  arising  from  the  scatter¬ 
ing  mechanism  and  the  antenna  patterns.  We  define  the  direction  to  the 
source  in  terms  of  an  angle,  a.  ,  viewed  from  the  receiver.  The 

weighted  Doppler  spectra,  P  (a,  ,f. )  ,  are  then  determined  by  applying 

w  r  i 

a  normalized  angular  scattering  function,  Wg(o^)  ,  and  a  combined  re¬ 
ceiver-transmitter  antenna  gain  function,  G  (ct^-a  )  ,  to  the  unweighted 
spectra,  I.e. 

PwO*r’fi^  =  P(ar,fi)  W^ou)  G^o^-a^)  (4.8) 


41 


Each  Doppler  frequency  is  converted  to  an  apparent  angular  source  direc¬ 
tion  through  use  of  the  Doppler  relation  In  (2.15),  which  for  a  uniform 


transverse  wind  velocity*,  vQ  ,  yields 


(4.9) 


On  the  basis  of  the  computed  angular  source  direction,  each  spectral  com¬ 
ponent  Is  weighted  in  accordance  with  the  classical  isotropic  turbulence 
model,  a  sin(x)/x  receiving  array  pattern,  and  a  Gaussian  transmitting 
pattern.  The  normalized  scattering  function  is  evaluated  for  a  given  ele¬ 
vation  angle,  3  »  using  the  following  expression  with  the  numerator  ob¬ 

tained  from  equation  2.20 


W3(oti) 


0.51  +  jVa-  +  e2 


-.-m 


0.84 


where:  m  =  11/3 
8  >  0.25° 


(4.10) 


The  combined  antenna  gain  for  the  assumed  patterns  is  evaluated  for  each 
source  direction  from  the  following 


Gc(V“r>  ■ 


2 


s1n[a(aj-a  )] 

air(a.-a  )/i80 

i  r 


exp[-  ba*] 


where:  a  =  549.900 
b  =  0.03429 


(4.11) 


The  sign  convention  for  the  wind  velocity  is  negative  for  a  SE  wind  and 
positive  for  a  NW  wind. 
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The  beam  constants  a  and  b  are  determined  for  half-power  beamwidths 
of  0.29°  and  3°,  respectively.  The  transmitting  beam  constant,  b  , 
includes  the  path  geometry  scaling  factor  which  permits  angular  measure¬ 
ments  to  be  referenced  to  the  receiver. 

4.3  Average  Doppler  Shift 

The  average  Doppler  shift  is  obtained  by  computing  the  first  moment 
(spectrum  centroid)  of  the  angular  Doppler  spectra.  For  the  unweighted 
spectra,  the  average  is  given  by 


100 

f  = 

i*l 

Similarly,  'or  the  weighted  spectra 
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(4.13) 


4.4  Average  Wind  Velocity 

The  average  transverse  wind  velocity  is  determined  using  the  results 
of  averaging  the  unweighted  and  weighted  Doppler  spectra.  An  unweighted 
average  velocity,  v  ,  is  obtained  directly  from  the  measured  unweighted 
average  Doppler  shift,  f  ,  by  using  equation  2.15,  i.e. 
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(4.14) 


A  weighted  average  wind  velocity,  v  ,  is  computed  to  include  correc- 

w 

tions  for  the  effects  of  the  scattering  mechanism  and  the  antenna  pat¬ 
terns.  This  result  is  obtained  from  a  linear  Interpolation  of  the  mea¬ 
sured  weighted  average  Doppler  shift,  f  ,  (for  uniform  wind  v  )  and 

W  Q 

the  measured  unweighted  average  Doppler  shift 


(4.15) 


In  these  expressions,  the  wind  speed  is  measured  in  meters  per  second 
and  all  angles  are  measured  in  degrees. 
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Chapter  5 


DISCUSSION  OF  DATA 


5.1  PIBAL  Mind  Data 

Local  wind  conditions  were  influenced  during  the  tests  by  a  low 
pressure  weather  system  centered  off  the  coast  of  the  state  of  Washing¬ 
ton.  About  six  hours  before  the  onset  (0100  April  21)  of  the  experiment, 
a  low  pressure  frontal  system  moved  through  the  area  with  slight  preci¬ 
pitation  and  some  gusty  winds.  During  the  two  11 -hour  periods  between 
0800  and  1800  when  pilot  balloons  (PIBALS)  were  launched,  the  measured 
cross-path  winds  were  generally  from  the  southeast.  Steady  southeast¬ 
erly  winds  prevailed  at  nearly  all  altitudes  on  the  first  day  of  the  ex¬ 
periment.  The  winds  were  moderate  and  progressively  increased  during 
the  day  reaching  maximum  velocities  in  mid-afternoon.  In  contrast, 
light  variable  winds  persisted  throughout  the  second  day  when  wind  speeds 
rarely  exceeded  two  meters  per  second.  Typical  transverse  winds  on  this 
day  were  less  than  one  meter  per  second  and  on  occasions  wind  direction 
changes  were  observed  between  SE  and  NW. 

PIBAL  readings  were  taken  every  30  seconds  at  approximately  100-m 
height  intervals.  PIBAL  data  above  the  radio  horizon  were  averaged  over 
400-m  vertical  segments  for  comparison  with  inferred  results  from  radio 
Doppler  measurements.  Successive  averages  were  taken  throughout  the 
vertical  extent  (1600  meters)  of  the  common  volume  region.  Average 
transverse  wind  velocities  in  the  four  vertical  segments  of  the  common 
volume  are  tabulated  in  Tables  5.1  and  5.2  for  the  PIBAL  wind  soundings 
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Table  5.1 


PIBAL  WIND-MEASUREMENT  DATA  ON  21  APRIL  1978* 


PIBAL 

Release 

Average  Transverse  Wind  Velocity  (m/s) 

1st  400-m  2nd  400-m  3rd  400-m  4th  400-m 

0805 

1.93 

SE 

3.58  SE 

5.69  SE 

7.12  SE 

0905 

1.06 

SE 

3.37  SE 

5.61  SE 

6.85  SE 

1005 

1.16 

SE 

2.73  SE 

5.19  SE 

5.86  SE 

1105 

1.84 

SE 

1.96  SE 

3.96  SE 

5.12  SE 

1205 

4.10 

SE 

4.49  SE 

1 .68  SE 

4.19  SE 

1305 

3.95 

SE 

4.18  SE 

4.55  SE 

4.02  SE 

1405 

3.21 

SE 

3.95  SE 

4.81  SE 

4.01  SE 

1505 

1.98 

SE 

1.08  SE 

3.51  SE 

4.51  SE 

1605 

3.52 

SE 

4.32  SE 

4.74  SE 

4.65  SE 

1705 

4.96 

SE 

4.65  SE 

2.51  SE 

3.55  SE 

1805 

3.15 

SE 

2.40  SE 

2.92  SE 

3.69  SE 

PIBAL 

PIBAL 

Release 

Table  5.2 

WIND-MEASUREMENT  DATA  ON  22  APRIL  1978* 

Average  Transverse  Wind  Velocity  (m/s) 

1st  400-m  2nd  400-m  3rd  400-m  4th  400-m 

0805 

0.10  NW 

0.92  NW 

0.48  SE 

2.82  SE 

0905 

0.12  NW 

1.49  NW 

0.14  SE 

1.76  SE 

1005 

0.34  SE 

0.09  SE 

0.53  NW 

0.98  SE 

1105 

0.54  NW 

0.14  NW 

1.00  NW 

0.08  SE 

1205 

0.43  NW 

0.78  SE 

0.99  NW 

1.13  NW 

1305 

0.09  SE 

0.29  SE 

2.17  NW 

0.40  NW 

1405 

0.53  SE 

1.88  SE 

0.56  NW 

2.25  NW 

1505 

1.45  SE 

0.73  SE 

1.07  NW 

0.09  NW 

1605 

1.01  SE 

0.72  SE 

0.98  NW 

0.52  NW 

1705 

0.45  SE 

0.22  SE 

0.70  NW 

2.72  NW 

1805 

0.89  NW 

0.18  NW 

2.17  NW 

1.02  NW 

* 

Average  transverse  wind  in 

400-m  vertical 

segments  above 

radio  horizoi 
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taken  respectively  on  April  21  and  22.  Examples  of  vertical  profiles  of 
the  cross-path  wind  component  observed  on  the  two  days  are  shown  in  Fig¬ 
ure  5.1  [(a)  thru  (f)].  Data  from  the  0805  launch  on  April  21  (Figure 
5.1a)  reflected  a  fairly  strong  wind  shear  between  the  1300  and  1600- 
meter  elevations.  Within  this  300-meter  vertical  range,  the  wind  changed 
from  0.8  m/s  NW  to  12  m/s  SE.  This  shear  dissipated  by  the  next  wind 
sounding  one  hour  later.  No  other  anomalous  wind  conditions  were  ob¬ 
served  for  the  duration  of  the  experiment. 

5.2  Wind-Produced  Doppler 

The  Doppler  spectra  presented  in  this  section  are  computed  from  600 
contiguous  sets  of  data  samples  representing  15  seconds  of  real-time  data. 
Each  data  set  is  comprised  of  twelve  amplitudes  and  twelve  relative 
phases.  Computations  are  made  from  the  same  data  sets  for  a  sequence  of 
array  pointing  angles  about  the  great-circle  bearing.  In  order  to  simul¬ 
taneously  measure  upwind  and  downwind  Doppler  shifts,  spectral  observa¬ 
tion  windows  are  centered  about  the  zero-Doppler  frequency.  Observation 
windows  used  in  the  spectrum  analysis  are  narrower  than  the  system  noise 
bandwidth  (23  Hz)  but  are  broad  enough  to  encompass  the  ambient  wind- 
produced  Doppler  shifts.  For  the  wind  conditions  experienced  during  the 
tests,  an  8  to  10-Hz  wide  observation  window  was  found  to  be  sufficient 
to  completely  envelop  the  wind-produced  Doppler  spectrum.  Outside  the 
selected  observation  windows,  the  Doppler  spectra  have  relatively  low 
power  spectral  densities.  Each  computed  spectrum  contains  100  Doppler 
frequencies  spaced  from  about  0.08  to  0.10  Hz  apart  depending  upon  the 
width  of  the  spectral  observation  window.  The  spectral  resolution 
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Fig.  5.1.  CONTINUED. 
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attainable  in  the  coherent  integration  process  of  the  600  data  samples 
is  0.07  Hz.  However,  the  realizable  resolution  is  more  on  the  order  of 
0.1  Hz  due  to  the  short-term  instabilities  (1.5  parts  in  10^)  of  the 
quartz  oscillators. 

5.2.1  Unweighted  Angular  Spectra 

Examples  of  some  unweighted  Doppler  spectra  are  shown  in  Figure  5.2 
[(a)  thru  (1)].  These  spectra  are  taken  for  symmetrical  (NW/SE)  beam 
pointing  directions  about  the  great-circle  bearing  at  0°  azimuth.  The 
spectra  are  plotted  over  a  40-dB  range  in  0.25°  increments  out  to  one 
degree  either  side  of  the  great-circle  bearing.  Two  15-second  samples 
taken  one  minute  apart  are  presented  for  each  hourly  example  in  order  to 
show  temporal  variations  in  the  wind-produced  Doppler  shifts. 

Initial  observations  of  the  high-resolution  spectra  show  the  pres¬ 
ence  of  a  coherent  stationary  source  with  a  very  narrow  spectral  width 
at  zero-Doppler  frequency.  This  spectral  component  is  most  likely  a  re¬ 
sult  of  knife-edge  diffraction  over  one  or  more  of  the  mountain  ridges 
between  the  transmitter  and  receiver.  It  is  seen  that  the  diffracted 
component  increases  in  intensity  as  the  array  scan  approaches  the  great- 
circle.  Occasionally,  when  the  array  beam  is  pointed  on-axis,  the  level 
of  the  diffracted  component  exceeds  the  40-dB  plotting  range  and  a  lin¬ 
ear  foldover  (overshoot)  occurs.  Aside  from  its  intensity  level,  the 
diffraction  component  is  totally  invariant  to  the  azimuthal  pointing 
direction  of  the  array. 

In  contrast,  the  wind-produced  Doppler  shifts  show  systematic  varia¬ 


tions  with  the  array  pointing  direction.  A  qualitative  analysis  of  the 
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Fig.  5.2.  CONTINUED. 


RELATIVE  POWER  (dB) 


Fig.  5.2.  CONTINUED. 


(h) 

Fig.  5.2.  CONTINUED. 


2.  CONTINUED. 


Fig.  5.2.  CONTINUED. 


FREQUENCY  (Hz) 


00 

Fig.  5.2.  CONTINUED. 


66 


wind-produced  spectra  shows  progressive  increases  in  the  average  Doppler 
shift  as  the  array  is  pointed  off-axis.  For  the  examples  shown  in  Figure 
5.2,  SE  winds  are  apparent  at  all  angles.  Doppler  shifts  are  generally 
positive  for  array  pointing  angles  in  the  SE  (upwind)  direction  and  nega¬ 
tive  for  angles  in  the  NW  (downwind)  direction.  At  most  off-axis  angles, 
the  upwind  and  downwind  Doppler  spectra  are  nearly  antisymmetric  about 
the  zero-Doppler  frequency  which  suggests  some  degree  of  spatial  uni¬ 
formity  in  the  transverse  winds.  On-axis  the  Doppler  spectra  are  gen¬ 
erally  symmetrical  about  zero  Doppler.  Any  asymmetries  appearing  in  the 
spectra  when  the  array  beam  is  aligned  along  the  great-circle  can  either 
be  attributed  to  a  vertical  wind  component  or  to  spatial  non-uniformities 
in  the  transverse  wind.  Birkemeier  et  al  (1977)  point  out  that  such 
spectral  asymmetries  also  result  from  the  misalignment  of  the  trans¬ 
mitting  and  receiving  antennas.  Misalignment  errors,  however,  are  elimi¬ 
nated  when  using  a  phased-array. 

Although  no  quantitative  analysis  is  presented,  other  observations 
from  Figure  5.2  reveal  definite  spectral  broadening  during  afternoon 
periods  when  the  atmosphere  is  fairly  well  mixed.  This  is  a  logical 
consequence  of  diffuse  (turbulent)  scattering  processes  which  produce 
wide  variations  in  the  signal's  angl e-of-arrival .  On  the  other  hand, 
in  the  morning  hours,  the  atmosphere  is  relatively  stable  and  is  more 
likely  to  be  partially  stratified.  In  such  situations,  a  quasi-coherent 
signal  is  received  that  arrives  via  a  number  of  specular  reflection 
points  distributed  over  a  small  angular  region  near  the  great-circle 
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bearing.  As  a  result  of  the  narrow  range  of  angles-of-arrival ,  the  sig¬ 
nal  spectra  are  not  nearly  as  broad  as  the  spectra  produced  by  signals 
propagating  through  a  turbulent  atmosphere.  The  spectra  are  also  seen 
to  be  somewhat  broader  on  the  first  day  [Figure  5.2  (a)  thru  (f)]  of  the 
experiment  than  on  the  second  day  [Figure  5.2  (g)  thru  (1)].  This  vari¬ 
ance  appears  directly  attributable  to  the  contrasting  wind  conditions  on 
the  two  days. 

The  short-term  variability  of  the  wind  can  be  observed  from  temp¬ 
oral  changes  in  the  Doppler  spectra.  Of  particular  importance  in  this 
regard  are  the  temporal  changes  in  the  average  Doppler  shift  as  a  func¬ 
tion  of  the  array  pointing  angle.  Figure  5.3  [(a)  thru  (f)]  shows  the 
measured  average  Doppler  shifts  for  the  spectra  given  in  Figure  5.2. 
Average  Doppler  shifts  are  computed  in  accordance  with  equation  4.12  for 
21  equally-spaced  (0.15°)  pointing  angles  between  1.5°  NW  and  1.5°  SE. 
The  slopes  of  these  curves  are  directly  proportional  to  the  prevailing 
wind  speeds.  Averages  taken  one  minute  apart  are  shown  on  the  same 
plot  for  comparative  analysis  of  short-term  variations.  The  15-second 
samples  are  generally  in  close  agreement  from  one  minute  to  the  next. 
This  suggests  reasonably  good  correlation  of  wind  conditions  within  one 
minute  intervals.  As  evidenced  by  the  near  linearity  of  the  Doppler 
shift  curves,  there  is  some  apparent  uniformity  in  the  wind  structure 
over  the  spatial  regions  probed.  Occasionally,  however,  there  are 
noticeable  departures  from  this  linear  trend  that  may  be  attributable 
to  variable  wind  conditions.  For  example,  in  Figure  5.3a  at  0803  on  21 
April  a  definite  change  in  the  polarity  of  the  average  Doppler  shift  is 
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Fig.  5.3.  CONTINUED. 
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observed  when  the  array  is  pointed  in  the  vicinity  of  0.3°  NW.  An  off- 
axis  polarity  change  such  as  this  is  reflective  of  a  corresponding  wind 
direction  change.  Supporting  this  are  the  wind-shear  conditions  re¬ 
flected  in  the  PIBAL  measurement  data  (Figure  5.1a)  at  the  same  hour. 
There  are  nonlinearities  in  the  average  Doppler  shift  curves  that  are 
not  necessarily  associated  with  the  wind.  These  nonlinearities  occur 
near  the  extremities  of  the  array  sector  scan  and  are  most  likely  due 
to  the  combined  off-axis  effects  of  atmospheric  scattering  and  secondary 
grating  lobes.  Effects  of  the  scattering  mechanism  can  be  taken  into 
account  by  weighting  the  Doppler  spectra  in  accordance  with  an  appropri¬ 
ate  atmospheric  scattering  model.  The  data  results  from  this  effort  are 
discussed  in  the  following  section. 

5.2.2  Weighted  Angular  Spectra 

The  Doppler  spectra  presented  in  the  previous  section  are  essen¬ 
tially  unweighted  in  that  the  spectral  components  are  derived  from  data 
samples  that  are  uniformly  weighted  over  the  angular  range  scanned  by 
the  array.  Here  we  take  into  account  the  apparent  source  direction  of 
each  spectral  component  and  the  angular  scattering  properties  of  the 
atmosphere.  Weighting  factors  based  on  angle-of-arrival  scattering 
variations  are  then  applied  to  the  individual  spectral  components  in  ac¬ 
cordance  with  procedures  outlined  in  section  4.2. 

Figure  5.4  [(a)  thru  (1)]  shows  the  results  of  applying  a  normal¬ 
ized  Tatarski  scattering  model  to  the  unweighted  spectral  data  in  Fig¬ 
ure  5.2.  These  spectra  are  obtained  for  an  assumed  uniform  wind  and  a 
turbulent  scattering  region  situated  at  an  elevation  angle  of  0.25°. 
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Fig.  5.4.  WEIGHTED  WIND-PRODUCED  DOPPLER  SPECTRA.  Angle-of-arrfval  spec¬ 
tra  are  weighted  In  accordance  with  (1)  a  normalized  Tatarskl  scattering 
model,  (2)  a  s1n(x)/x  receiving  array  pattern,  and  (3)  a  Gaussian  trans¬ 
mitting  pattern  under  the  assumption  of  a  uniform  wind  and  a  source  ele¬ 
vation  angle,  8  *  0.25°. 
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